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Preface 


The Compag Portable Mathematics Library (CPML) is a set of mathematical 
routines that are accessed from high-level languages (such as Fortran and C) 
which support mathematical functions. Many CPML routines can also be called 
directly using standard call interfaces, but it is recommended that you invoke 
CPML routines only from a high-level language. 


Intended Audience 
This book is for compiler writers, system programmers, and application 
programmers who want to use CPML routines. 

Document Structure 


This manual consists of the following: 


Chapter 1 gives a general overview of the mathematics library and discusses 
supported data types, exception behavior, and IEEE considerations. 


Chapter 2 explains the presentation format of a CPML routine and how to 
interpret a routine’s interface. It also alphabetically lists and describes the 
routines. 


Appendix A lists the floating-point boundary values used by the CPML routines. 
Appendix B contains the complete list of entry-point names. 


The Glossary lists mathematical terms and symbolic names used in this manual, 
and provides a brief definition. 


Related Documents 


Some books in Compaq’s documentation sets help meet the needs of several 
audiences. For example, the information in some system books is also used by 
programmers. Keep this in mind when searching for information on specific 
topics. 


Use the documentation overview and the master index information for your 
operating system when searching for hardcopy information on a topic. They 
provide information on all of the books in your operating system's documentation 
set. 

CPML Documentation 

For additional information about CPML, you can access the Compaq CPML 
website at the following location: 


http: //www.compagq.com/math 


vii 


OpenVMS Documentation 


For additional information about Compag OpenVMS products and services, access 
the Compaq website at the following location: 


http: //www.openvms.compaq.com/ 


How to Order Additional Documentation 


Use the following World Wide Web address to order additional documentation: 
http: //www.openvms.compaq.com/ 


If you need help deciding which documentation best meets your needs, call 
800-282-6672. 


Reader’s Comments 


Compaq welcomes your comments about this manual. Please send comments to 
one of the following: 


Internet openvmsdoc@compaq.com 


Mail Compaq Computer Corporation 
OSSG Documentation Group, ZK O3-4/U 08 
110 Spit Brook Rd. 
Nashua, NH 03062-2698 


Conventions 


viii 


In this book, every use of OpenVMS means Compaq’s OpenVMS operating 
system, and every use of UNIX means Compaq’s Tru64 UNIX operating system. 


The following conventions are used in this manual: 


Ctrl/x A sequence such as Ctrl/x indicates that you must hold down 
the key labeled Ctrl while you press another key or a pointing 
device button. 


PF1x A sequence such as PF 1 x indicates that you must first press 
and release the key labeled PF 1 and then press and release 
another key or a pointing device button. 


Return In examples, a key name enclosed in a box indicates that 
you press a key on the keyboard. (In text, a key name is not 
endiosed in a box.) 


In the HTML version of this document, this convention appears 
in brackets, rather than a box. 


A horizontal ellipsis in examples indicates one of the following 

possibilities: 

¢ Additional optional arguments in a statement have been 
omitted. 


¢ The preceding item or items can be repeated one or more 
times. 


¢ Additional parameters, values, or other information can be 
entered. 


A vertical ellipsis indicates the omission of items from a code 
example or command format; the items are omitted because 
they are not important to the topic being discussed. 


[] 


Ace 


bold text 


italic text 


UPPERCASE TEXT 


Monospace text 


numbers 


In command format descriptions, parentheses indicate that you 
must enclose choices in parentheses if you specify more than 
one. 


In command format descriptions, brackets indicate optional 
choices. You can choose one or more items or no items. 

Do not type the brackets on the command line. However, 
you must include the brackets in the syntax for OpenVMS 
directory specifications and for a substring specification in an 
assignment statement. 


In command format descriptions, vertical bars separate choices 
within brackets or braces. Within brackets, the choices are 
optional; within braces, at least one choice is required. Do not 
type the vertical bars on the command line. 


In command format descriptions, braces indicate required 
choices; you must choose at least one of the items listed. Do 
not type the braces on the command line. 


This typeface represents the introduction of a new term. It 
also represents the name of an argument, an attribute, or a 
reason. 


Italic text indicates important information, complete titles 
of manuals, or variables. Variables include information that 
varies in system output (Internal error number), in command 
lines ((PRODUCER=name), and in command parameters in 
text (where dd represents the predefined code for the device 
type). 


Uppercase text indicates a command, the name of a routine, 
the name of a file, or the abbreviation for a system privilege. 


Monospace type indicates code examples and interactive 
screen displays. In the C programming language, monospace 
type in text identifies the following elements: keywords, the 
names of independently compiled external functions and files, 
syntax summaries, and references to variables or identifiers 
introduced in an example. 


A hyphen at the end of a command format description, 
command line, or code line indicates that the command or 
statement continues on the following line 


All numbers in text are assumed to be decimal unless 
otherwise noted. Nondecimal radixes—binary, octal, or 
hexadecimal—are explicitly indicated. 


1 


Introduction to CPML 


The Compaq Portable Mathematics Library (referred to as CPML) includes a 
wide variety of mathematical routines that cover the following areas: 


¢ Floating-point trigonometric function evaluation 

¢ Exponentiation, logarithmic, power function evaluation 
¢ Hyperbolic function evaluation 

e Algebraic function evaluation 

* Complex function evaluation 

¢« Complex exponentiation 

e Miscellaneous function evaluation 


This manual documents the CPML routines and, in particular, how they behave 
when given an exceptional input argument. It also documents operating system 
entry points and supported floating-point data types. 


1.1 Overview 


Developing software within the confines of high-level languages like Fortran 

and C greatly increases the portability and maintainability of your source code. 
Many high-level languages support mathematical function evaluation. CPML was 
developed to provide a common set of routines that supports many of the common 
mathematical functions across a wide variety of operating systems, hardware 
architectures, and languages. 


In most cases, the common mathematical functions behave in the same way 

for all languages and platforms. Occasionally, however, high-level language 
definitions of the same mathematical function differ for specific input values. For 
example, in Fortran, log(-1.0) causes a program abort, while in C, log(-1.0) quietly 
returns a system-defined value. 


This document uses the term exceptional arguments to refer to values in the 
following situations: 


e Values for which high-level languages disagree on the function behavior 
¢« Values that are mathematically undefined or out of range 

e Values for which the function would overflow or underflow 

See Section 1.3 for more detail on exceptional arguments. 


To provide uniform quality of mathematical functions for all languages on your 
system, CPML traps exceptional arguments and invokes a system-specific routine 
called the CPML exception handler. The exception handler is designed to work 
with high-level language compilers and run-time libraries (RTLS) to provide 
specific language semantics for exceptional arguments. This means that the 
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user-visible behavior of a given function called from a given language is not 
necessarily determined by the routines in the CPML library but rather by a 
combination of several entities acting in concert. 


Note 


Compaq strongly recommends that you limit your access to the CPML 
routines documented in this manual to the high-level language syntax 
of your choice, thereby guaranteeing the behavior of the routines across 
platforms. Because of the complex relationship between high-level 
languages and CPML routines, the behavior of direct calls to CPML 
routines may change from release to release. 


1.2 Data Types 
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CPML is designed to support mathematics function evaluation for multiple data 
types. These data types include integer, floating-point, and complex floating-point. 


The integer data type, identified as int throughout this manual, is the natural 
size signed integer for a particular platform. On a 32-bit system, int is a 32-bit 
signed integer, and on a 64-bit system, int is a 64-bit signed integer. 


The floating-point types referred to in this document are F_FLOAT, G FLOAT, 
X_FLOAT, S_FLOAT, and T_FLOAT, respectively. When it is not necessary to 
distinguish between the different floating types, they are referred to collectively 
as F_TYPE. Your platform may support all or a subset of these floating-point data 
types. For example, CPML on OpenVMS Alpha systems supports the following 
floating-point data types: VAX single: and double-precision, IEEE single and 
double-precision, and IEEE extended-precision. CPML on Compaq Tru64 UNIX 
Alpha systems supports only IEEE single: and double-precision data types. 
Table 1-1 describes the floating-point data types. 


Table 1-1 Floating-Point Data Types 


F_TYPE Description 

S_ FLOAT 32-bit IEEE single-precision number 
T_FLOAT 64-bit IEEE double-precision number 
X_FLOAT 128-bit IEEE extended-precision number 
F_FLOAT 32-bit VAX single-precision number 

G_ FLOAT 64-bit VAX double-precision number 


In addition to the data types mentioned in Table 1-1, CPML also provides 
routines that return two values of the same floating-point type, for example, two 
S_TYPE values or twoG TYPE values. In the discussion that follows, these 
pairs of floating-point data type values are referred to as F COMPLEX. Refer 
to Table 1-2. This document uses F_COMPLE-X to indicate that a given routine 
returns two different values of the same floating-point data type. 


The mechanism for returning two floating-point values from CPML routines 
varies from platform to platform. However, on OpenVMS Alpha systems, 
F_ COMPLEX data is returned in consecutive floating-point registers and is 
accessible only through a high-level language, like Fortran, that specifically 
allows access to it. 
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A complex number, z, is defined as an ordered pair of real numbers. The 
convention used in this manual to define an ordered pair of real numbers as 
complex is as follows: 


e The first number is the real part of the complex number. 


« The second number is preceded by i and is the imaginary part of the complex 
number. 


e« A separator character (plus sign) is used to associate and separate the real 
and the imaginary number. 


For example: 


Z =x +iy 
z =sin xX +icos y 


CPML includes complex functions, for example, the complex sine, csin(x,y), 
defined to be sin(x +iy). Complex function routines like csin(), which have 
complex input, accept floating-point numbers in pairs and treat them as if they 
are real and imaginary parts of a complex number. 


In the previous two examples, the first floating-point values are defined by 

xX and sin x, respectively, and are the real part of the complex number. The 
second floating-point values used in the examples are defined by iy and icos y, 
respectively, and are the imaginary part of the complex number. Similarly, CPML 
routines that return complex function values return two floating-point values. 
Taken together, these two floating-point values represent a complex number. 


CPML supports the floating-point complex types described in Table 1-2. You can 
access CPML complex functions only through high-level languages that support 
the complex data type. Use only the data types supported by your system. 


Table 1-2 Floating-Point Complex Data Types 
F_COMPLEX Description’ 


S_FLOAT_COMPLEX An ordered pair of S FLOAT quantities, representing a single 
precision complex number 


T_FLOAT_COMPLEX An ordered pair of T FLOAT quantities, representing a double 
precision complex number 


X_FLOAT_COMPLEX An ordered pair of X_FLOAT quantities, representing an 
extended-precision complex number 


F_FLOAT_COMPLEX An ordered pair of F_FLOAT quantities, representing a single 
precision complex number 


G_FLOAT_COMPLEX An ordered pair of G FLOAT quantities, representing a double- 
precision complex number 


1The lower addressed quantity is the real part; the higher addressed quantity is the imaginary part. 


1.3 Exceptional Arguments 


Not all mathematical functions are capable of returning a meaningful result 

for all input argument values. Any argument value passed to a CPML routine 
that does not return a meaningful result, or is defined differently for different 
environments, is referred to as an exceptional argument. Exceptional arguments 
that result in an exception behavior are documented in the Exceptions section of 
each CPML routine in Chapter 2. 
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Exceptional arguments typically fall into one of two categories: 


¢ Domain errors or invalid arguments. These are arguments for which a 
function is not defined. For example, the inverse sine function, asin, is 
defined only for arguments between -1 and +1 indusive. Attempting to 
evaluate acos(-2) or acos(3) results in a domain error or invalid argument 
error. 


e Range errors. These errors occur when a mathematically valid argument 
results in a function value that exceeds the range of representable values for 
the floating-point data type. Appendix A gives the approximate minimum and 
maximum values representable for each floating-point data type. 


1.4 Exception Conditions and Exception Behavior 
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CPML routines are designed to provide predictable and platform-consistent 
exception conditions and behavior. When an exception is triggered in a CPML 
routine, two pieces of information can be generated and made available to the 
calling program for exception handling: 


¢ A notification that an exception has occurred. The mechanics of exception 
notification vary from platform to platform (for example, signaling, trapping, 
set errno). 


e A return value. If your environment allows your routine to continue after 
raising an exception condition (with an exception handler for example), then 
a return value is made available upon completion of the routine. 


The exception condition-handling mechanisms on your platform dictate how you 
can recover from an exception condition, and whether you can expect to receive 
an exception notification, a return value, or both, from a CPML routine. 


The Exceptions section of each CPML routine documents each exceptional 
argument that results in an exception behavior. In addition to the exceptional 
arguments, an indication of how the CPML routines treat each argument is given. 
Exceptional arguments are sometimes presented in terms of symbolic constants. 


For example, the following table lists the exceptional arguments of the 
exponential routine, exp(x): 


Exceptional Argument Exception Condition/Routine Behavior 


x >In(max_float) Overflow 
x <In(min_float) Underflow 


The exceptional arguments indicate that whenever x >In(max_float) 
or x <In(min_float), CPML recognizes an overflow or underflow condition, 
respectively. 


The symbolic constants In(max_float) and In(min_float) represent the natural 
log of the maximum and minimum representable values of the floating-point 
data type in question. The actual values of In(max_float) and In(min_float) are 
described in Appendix A. 


CPML recognizes three predefined conditions: overflow, underflow, and invalid 
argument. Table 1-3 describes the default action and return value of each 
condition. 
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Table 1-3 Default Action and Return Values for Exception Conditions 


Exception Condition Default Action Return Value 
Overflow Trap HUGE_RESULT 
U nderflow Continue Quietly ) 

Invalid argument Trap INV_RESULT 


The values HUGE RESULT and INV_RESULT are data-type dependent. 


For IEEE data types, HUGE RESULT and !INV_RESULT are the floating-point 
encodings for Infinity and NaN, respectively. 


For VAX data types, HUGE RESULT and INV_RESULT are max_float and 0, 
respectively. 


1.5 IEEE Std 754 Considerations 


The Institute of Electrical and Electronics Engineers (IEEE) ANSIA,EEE Std 
754-1985, IEEE Standard for Binary Floating-Point Arithmetic data types 
include denormalized numbers (very close to zero). The standard supports the 
concept of “Not-a-Number” or NaN to represent indeterminate quantities, and 
uses plus infinity or minus infinity (so that they behave in arithmetic) like the 
mathematical infinities. Whenever a CPML routine produces an overflow or 
indeterminate condition, it generates an infinity or NaN value. 


All CPML routines, except one, return a NaN result when presented with a NaN 
input. The only exception is pow(NaN,0) =1in ANSI C. 


1.6 X/Open Portability Guide Considerations 


Table 1-4 lists the routines described in this manual that conform to the 
requirements of the X/ Open Portability Guide, Version 4 (XPG4), or are 
implemented as UNIX extensions to the XPG4 standard (XPG4-UNIX). 
Descriptions of these routines appear in Chapter 2 under the generic 
function name listed in Table 1-4. Platform-specific entry-points are listed 
in Appendix B. 


Table 1-4 XPG4 Conformant Routines 


Routine Conforms to Standard Generic Function Name 
acos XPG4 acos 

acosh XPG4-UNIX acosh 

asin XPG4 asin 

asinh XPG4-UNIX asinh 

atan XPG4 atan 

atan2 XPG4 atan 

atanh XPG4-UNIX atanh 

ceil XPG4 ceil 

cos XPG4 cos 


(continued on next page) 
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Table 1—4 (Cont.) XPG4 Conformant Routines 


Routine Conforms to Standard Generic Function Name 
cosh XPG4 cosh 

cot XPG4 cot 

erf XPG4 erf 

erfc XPG4 erf 

exp XPG4 exp 
expm1 XPG4-UNIX exp 
fabs XPG4 fabs 
floor XPG4 floor 
fmod XPG4 fmod 
frexp XPG4 frexp 
gamma XPG4 Igamma 
hypot XPG4 hypot 
ilogb XPG4-UNIX ilogb 
isnan XPG4 isnan 

jo XPG4 bessel 
jl XPG4 bessel 
jn XPG4 bessel 
Idexp XPG4 Idexp 
Igamma XPG4 Igamma 
log XPG4 log 
log10 XPG4 log 
loglp XPG4-UNIX log 
logb XPG4-UNIX logb 
modf XPG4 modf 
nextafter XPG4-UNIX nextafter 
pow XPG4 pow 
remainder XPG4-UNIX remainder 
rint XPG4-UNIX rint 
scalb XPG4-UNIX scalb 
sin XPG4 sin 

sinh XPG4 sinh 
tan XPG4 tan 
tanh XPG4 tanh 

yO XPG4 bessel 
yl XPG4 bessel 
yn XPG4 bessel 
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CPML Routines 


CPML routines can be accessed from high-level languages that support 
mathematical functions (such as Fortran and C), or called directly using standard 
call interfaces. It is highly recommended that you invoke CPML routines only 
from a high-level language. 


CPML routines are documented with generic names, and with the symbol 
F_TYPE to indicate generic floating-point values (e.g. F_ TYPE sqrt (F_TYPE x)). 


To determine the appropriate names and interfaces within a specific programming 
language (e.g. float sqrtf(float x) or REAL*4 SQRT), refer to that language's 
documentation. 


To enable the use of CPML routines which are not provided by your high-level 
language, the actual CPML entrynames are provided. 


Note: CPML routines which return complex numbers ("F_ COMPLEX") usea 
private interface. Therefore, they can only be called from high-level languages 
that support that interface. 


The Data Types S FLOAT, T FLOAT and X_FLOAT refer to |EEE format 
floating-point numbers of single-, double-, and quad-precision, respectivel y. 
F_FLOAT and G FLOAT refer to VAX format single precision, and G-floating 
double-precision floating point numbers, respectively. 


For each CPML routine, "exceptional" input values are also provided. That is, 
values for which the function is mathematically undefined, or for which the 
output would be out of range for the floating-point type. 


Refer to your language’s documentation for information about how exceptions 
manifest themselves and how to control exception behavior. 


Further information is also available at the Compaq Math website at: 
http://www.compaq.com/math. 
2.1 CPML Routine Descriptions 


CPML routines are described in detail at the end of this chapter. Each CPML 
routine documented in this chapter is presented in the following format: 


e« Routine name—A brief name to identify the function of the routine. A routine 
may contain more than one function. 


e |nterface—What the routine expects to receive and what it returns. See 
Section 2.2 for more information. 


¢ Description—Additional information, including the permitted range of input 
values and generic calculations used to compute the results. 


e Exceptions—A description of how the routine behaves when given a specific 
exceptional input argument. 
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2.2 CPML Routine Interface 


The interface to each function is: 
RETURN_TYPE generic _interface_name (INPUT_ARG TYPE...) 
Each of these is described below. 


RETURN_TYPE 

The data type of the value that the routine returns to your application program. 
Each routine returns a specific class of data type. For example, either F_TYPE 
or F_ COMPLEX can appear in a CPML interface as described in Chapter 2. The 
supported data types are described in Section 1.2. 


generic_interface_name 

The generic name. CPML routines in this chapter are listed in alphabetic order 
by their interface names. Some CPML routines may be available in the syntax 
of your high-level language. Fortran and C are examples. To maximize the 
portability of your application, use the corresponding mathematical routine 
described in your high-level language, and directly call only the routines 
documented in this manual that are not supported by your language. Refer to 
Appendix B for the specific entry-point names needed to directly call a CPML 
routine from your platform. 


INPUT_ARG_TYPE... 

The number and type of input arguments provided by your application. Some 
routines require more than one argument. Arguments must be coded in the order 
shown in the interface section of each routine described in this chapter. The 
supported data types for arguments are described in Section 1.2. 


Note 


Unless otherwise noted, arguments are read-only and passed by value. 
Arguments passed by another mechanism are prefaced by an asterisk (*); 
for example, *n in the frexp() routine. 


2.3 Specific Entry-Point Names 


Each generic interface name documented in the interface section of a routine 
description corresponds to one or more specific entry-point names described 

in Appendix B. For example, on OpenVMS Alpha systems, the acosd function 
has five entry-point names, one for each available floating-point data type. 

The acosd entry-point names are math$acosd_f, math$acosd_s, math$acosd_x, 
math$acosd_g, and math$acosd_t. On Compaq Tru64 UNIX Alpha systems, 
the acosd function has two entry-point names corresponding to their supported 
data types: S FLOAT and T_FLOAT. The two entry-point names are acosdf for 
S FLOAT input arguments and acosd for T FLOAT arguments. Use the specific 
entry-point name that corresponds to the input argument data type. 
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2.4 Working with Exception Conditions 


Each CPML routine description contains a table of exceptions. Each exception 
listed in the table represents an exceptional case that is handled in a platform- 
specific manner. For example, the atan2() exception table contains the following 


two entries: 

Exceptional Argument Routine Behavior 
y=x=0 Invalid argument 
| y| =| x| =infinity Invalid argument 


The first entry describes an exception condition containing two input arguments 
with zero values. Upon detecting this error, the routine behavior signals the 
“invalid argument” condition. The second entry is applicable only to platforms 
supporting signed or unsigned infinity values. Here, if the absolute value of both 
input arguments is equal to infinity, an “invalid argument” condition is signaled. 


The exact behavior of a routine that detects an exceptional argument varies from 
platform to platform and is sometimes dependent on the environment in which 

it is called. The behavior you see depends on the platform and language used. 

It also depends on how the routine was called and the interaction of the various 
layers of software through which the call to the routine was made. Remember, 
access to a CPML routine can be made either through direct access (a CALL 
statement written by a programmer in a source code statement) or through 
indirect access (from compiler-implemented mathematical syntax). 


The default behavior for detecting the x=y=0 arguments is to generate an 
exception trap when accessing atan2() indirectly through Fortran compiler 
syntax. C compiler syntax for the atan2() routine sets errno and returns a NaN 
when give the same input. In these cases, your compiler documentation provides 
you with information on how to work with exception conditions. 


2.5 CPML Routine Interface Examples 


This section discusses the atan2() and cdiv() interfaces and explains how to 
interpret them. The explanations given in this section apply to all CPML 
routines. 


2.5.1 atan2() Interface 
The interface to the atan2() routine is: 


F_TYPE atan2 (F_TYPE y, F_TYPE x) 


The routine name atan2() is the high-level language source-level name that gets 
mapped to a specific entry-point name documented in Appendix B. This is the 
name that appears in compiler documentation for this mathematical routine. The 
appropriate entry-point name is automatically selected when atan2() is called 
from high-level language syntax. This selection depends upon the data type of 
the input arguments. If you make direct calls to this routine, you must manually 
select the proper entry-point name documented in Appendix B for the data type 
of your input arguments. 


The format of the atan2() routine shows that it expects to receive two input 
arguments by value. Both arguments must be the same F_TYPE. The returned 
value will also be the same F_TYPE as the input arguments. 
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For example, on OpenVMS Alpha systems, the G FLOAT entry-point name is 
math$atan2_g(). It takes two G FLOAT arguments by value and returns a 
G_ FLOAT result. 


For Compaq Tru64 UNIX Alpha systems, the S FLOAT entry-point name is 
atan2f(). The routine takes two S FLOAT input arguments by value and returns 
an S FLOAT result. 


2.5.2 cdiv() Interface 
The interface to the cdiv() routine is: 
F_ COMPLEX cdiv (F_TYPE a, F_TYPE b, F_TYPE c, F_TYPE d) 


The routine name cdiv() is the generic name that gets mapped to a specific entry- 
point name documented in Appendix B. Selection of the appropriate entry-point 
name is done automatically when cdiv() is called from high-level language syntax. 
This selection depends upon the data type of the input arguments. Again, if you 
make direct calls to this routine, you must manually select the proper entry-point 
name documented in Appendix B for the data type of your input arguments. 


The format of the cdiv() routine shows that it expects to receive four input 
arguments by value. All arguments must be the same F_TYPE. The returned 
value will be an F_COMPLEX data type and will be the same base data type as 
the input arguments. 


For example, on OpenVMS Alpha systems, the F_FLOAT entry-point name is 
math$cdiv_f(). This routine takes four F_FLOAT input arguments by value 
and returns an F_FLOAT_ COMPLEX result in an ordered pair of F_FLOAT 
quantities. 


For Compaq Tru64 UNIX Alpha systems, the S FLOAT entry-point name is 
cdivf(). This routine takes four S FLOAT input arguments by value and returns 
an S FLOAT_COMPLEX result. 
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acos - Arc Cosine of Angle 


Interface 
F_TYPE acos (F_TYPE x) 
F_TYPE acosd (F_TYPE x) 

Description 
acos() computes the principal value of the arc cosine of x in the interval [0,pi] 
radians for x in the interval [-1,1]. 
acosd() computes the principal value of the arc cosine of x in the interval [0,180] 
degrees for x in the interval [-1,1]. 

Exceptions 


Exceptional Argument Routine Behavior 


| x] >1 Invalid argument 
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acosh - Hyperbolic Arc Cosine of Angle 


Interface 
F_TYPE acosh (F_TYPE x) 

Description 
acosh() returns the hyperbolic arc cosine of x for x in the interval [1,+infinity]. 
acosh(x) =In(x +sqrt(x**2 - 1)). 
acosh() is the inverse function of cosh(). The definition of the acosh() function is 
acosh(cosh (x)) =x. 

Exceptions 


Exceptional Argument Routine Behavior 


x<l Invalid argument 
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asin - Arc Sine of Angle 


Interface 
F_TYPE asin (F_TYPE x) 
F_TYPE asind (F_TYPE x) 

Description 
asin() computes the principal value of the arc sine of x in the interval [-pi/2,pi/2] 
radians for x in the interval [-1,1]. 
asind() computes the principal value of the arc sine of x in the interval [-90,90] 
degrees for x in the interval [-1,1]. 

Exceptions 


Exceptional Argument Routine Behavior 


| x] >1 Invalid argument 
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asinh - Hyperbolic Arc Sine of Angle 


Interface 
F_TYPE asinh (F_TYPE x) 
Description 
asinh() returns the hyperbolic arc sine of x for x in the interval 
[-infinity, Hnfinity]. asinh(x) =In(x +sqrt(x**2 +1)). 
asinh() is the inverse function of sinh(). asinh(sinh (x)) =x. 
Exceptions 


None. 
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atan - Arc Tangent of Angle with One Argument 


Interface 
F_TYPE atan (F_TYPE x) 
F_TYPE atand (F_TYPE x) 

Description 
atan() computes the principal value of the arc tangent of x in the interval 
[-pi/2,pi/2] radians for x in the interval [-infinity, tinfinity]. 
atand() computes the principal value of the arc tangent of x in the interval 
[-90,90] degrees for x in the interval [-infinity, +infinity]. 

Exceptions 


None. 
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atan2 - Arc Tangent of Angle with Two Arguments 


Interface 


Description 


Exceptions 


F TYPE atan2 (F_TYPE y, F_TYPE x) 
F_TYPE atand2 (F_TYPE y, F_TYPE x) 


atan2() computes the angle in the interval [-pi,pi] whose arc tangent is y/x 
radians for x and y in the interval [-infinity, tinfinity]. The sign of atan2() is the 
same as the sign of y. The atan2(y, x) function is computed as follows, where f is 
the number of fraction bits associated with the data type: 


Value of Input Arguments Angle Returned 
x=0o0or gee n/2% (St gry) 
z>oOand y/x < gf atan (y/=) 

z<oOand ,/x < gftl T * (St gry) + atan (y/=) 


atand2() computes the angle in the interval [-180,180] whose arc tangent is y/x 
degrees for x and y in the interval [-infinity, +infinity]. The sign of atand2() is 
the same as the sign of y. 


Exceptional Argument Routine Behavior 
y=x=0 Invalid argument 
| y| =infinity and | x| =infinity Invalid argument 
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atanh - Hyperbolic Arc Tangent of Angle 


Interface 
F_TYPE atanh (F_TYPE x) 


Description 


atanh() returns the hyperbolic arc tangent of x for x in the interval (-1,1). atanh() 
is the inverse function of tanh(). atanh(tanh (x)) =x. 


atanh(x) is computed as 1/2 In((1+x)/(1-x)). 


Exceptions 


Exceptional Argument Routine Behavior 


| x| >or=1 Invalid argument 
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bessel - Bessel Functions 


Interface 
F_TYPE jO (F_TYPE x) 
F_TYPE j1 (F_TYPE x) 
F_TYPE jn (int n, F_TYPE x) 
F_TYPE yO (F_TYPE x) 
F_TYPE y1 (F_TYPE x) 
F_TYPE yn (int n, F_TYPE x) 


Description 


jO() and j1() return the value of the Bessel function of the first kind of orders 0 
and 1, respectively. 


jn() returns the value of the Bessel function of the first kind of order n. 


yO() and y1() return the value of the Bessel function of the second kind of orders 
0 and 1, respectively. 


yn() returns the value of the Bessel function of the second kind of order n. 
The value of x must be positive for the y family of Bessel functions. The value of 
n specifies some integer value. 


Exceptions 


Exceptional Argument Routine Behavior 
(yO(), y1(), yn()) x <O Invalid argument 
(yO(), y1(), yn()) x =0 Overflow 


The j1() and jn() functions can result in an underflow as x becomes small. The 
largest value of x for which this occurs is a function of n. 


The y1() and yn() functions can result in an overflow as x becomes small. The 
largest value of x for which this occurs is a function of n. 
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cabs - Complex Absolute Value 


Interface 
F_TYPE cabs (F_TYPE x, F_TYPE y) 
Description 
cabs(x,y) is defined as the square root of (x**2 +y**2) and returns the same value 
as hypot(x,y). 
Exceptions 
Exceptional Argument Routine Behavior 
sqrt(x**2 + y**2) >max_float Overflow 
See Also 


Appendix A, Critical Floating-Point Values 
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cbrt - Cube Root 


Interface 
F_TYPE cbrt (F_TYPE x) 


Description 


cort() returns the cube root of x. 


Exceptions 


None. 
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ccos - Cosine of Angle of a Complex Number 


Interface 
F_COMPLEX ccos (F_TYPE x, F_TYPE y) 
Description 
ccos() returns the cosine of a complex number, x +iy. 
ccos(x,y) is defined as cos (x +iy) =(cos x * cosh y - i * sin x * sinh y). 
Exceptions 
Exceptional Argument Routine Behavior 
| x| =infinity Invalid argument 
(sin x sinh y) > max_float Overflow 
(cos x cosh y) >max_float Overflow 
See Also 


Appendix A, Critical Floating-Point Values 
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cdiv - Complex Division 


Interface 
F_COMPLEX cdiv (F_TYPE a, F_TYPE b, F_TYPE c, F_TYPE d) 


Description 


cdiv() returns the quotient of two complex numbers: (a +ib)/(c +id). 


Exceptions 


Exceptional Argument Routine Behavior 


c=d=0 Invalid argument 


The quotient may overflow. 
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ceil - Ceiling 


Interface 
F_TYPE ceil (F_TYPE x) 


Description 


ceil() returns the smallest floating-point number of integral value greater than or 
equal to x. 


Exceptions 


None. 
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cexp - Complex Exponential 


Interface 
F_COMPLEX cexp (F_TYPE x, F_TYPE y) 
Description 
cexp() returns the exponential of a complex number. 
cexp(x,y) is defined as e**(x +iy) =e**x cos y +ie**x sin y. 
Exceptions 
Exceptional Argument Routine Behavior 
| y| =infinity Invalid argument 
| e**x cos y|_ >max_float Overflow 
| ex sin y|  >max_float Overflow 
See Also 


Appendix A, Critical Floating-Point Values 
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clog - Complex Natural Logarithm 


Interface 
F_COMPLEX clog (F_TYPE x, F_TYPE y) 


Description 
clog() returns the natural logarithm of a complex number. 
clog(x,y) is defined as In(x +iy) =1/2 In(x**2 + y**2) +i * atan2(y,x). 


Exceptions 


Exceptional Argument Routine Behavior 
y=x=0 Invalid argument 
| y| =| x| Snfinity Invalid argument 
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cmul - Complex Multiplication 


Interface 
F_COMPLEX cmul (F_TYPE a, F_TYPE b, F_TYPE c, F_TYPE d) 


Description 
cmul() returns the product of two complex numbers. 
cmul(a,b,c,d) is defined as (a +ib) * (c +id). 


Exceptions 


None. 


CPML-20 CPML Routines 


CPML Routines 
copysign - Copy Sign 


copysign - Copy Sign 


Interface 
F_TYPE copysign (F_TYPE x, F_TYPE y) 


Description 


copysign() returns x with the same sign as y. IEEE Std 754 requires 
copysign(x,NaN) = +x or -x. 


Exceptions 
None. 
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cos - Cosine of Angle 


Interface 
F_TYPE cos (F_TYPE x) 
F_TYPE cosd (F_TYPE x) 
Description 
cos() computes the cosine of x, measured in radians. 
cosd() computes the cosine of x, measured in degrees. 
Exceptions 


Exceptional Argument Routine Behavior 


| x] =infinity Invalid argument 


CPML-22 CPML Routines 


CPML Routines 
cosh - Hyperbolic Cosine of Angle 


cosh - Hyperbolic Cosine of Angle 


Interface 
F_TYPE cosh (F_TYPE x) 


Description 
cosh() computes the hyperbolic cosine of x. 
cosh(x) is defined as (exp(x) + exp(-x))/2. 


Exceptions 


Exceptional Argument Routine Behavior 


| x] >In(2 * max_float) Overflow 


See Also 
Appendix A, Critical Floating-Point Values 
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cot - Cotangent of Angle 


Interface 
F_TYPE cot (F_TYPE x) 
F_TYPE cotd (F_TYPE x) 
Description 
cot() computes the cotangent of x, measured in radians. 
cotd() computes the cotangent of x, measured in degrees. 
Exceptions 


Exceptional Argument Routine Behavior 
(cot) x=0 Overflow 
(cotd) | x| =multiples of 180 degrees Overflow 
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cpow - Complex Power 


Interface 
F_COMPLEX cpow (F_TYPE a, F_TYPE b, F_TYPE c, F_TYPE d) 
Description 
cpow() raises a complex base (a +ib) to a complex exponent (c +id). 
cpow(a,b,c,d) is defined as e**((c +id) In(a +ib)). 
Exceptions 
Exceptional Argument Routine Behavior 
sqrt (a**2 +b**2) >max_float Overflow 
c/2 * In(a**2 + b**2) >max_float Overflow 
c/2 * In(a*¥*2 + b**2) - (d * atan2(b,c)) >max_float Overflow 
a=b=c=d=0 Invalid argument 
See Also 


Appendix A, Critical Floating-Point Values 
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csin - Sine of Angle of a Complex Number 


Interface 
F_COMPLEX csin (F_TYPE x, F_TYPE y) 
Description 
csin( ) computes the sine of a complex number, x + iy. 
csin(x,y) is defined as csin (x +iy) =sin x * cosh y +i * cos x * sinh y. 
Exceptions 
Exceptional Argument Routine Behavior 
| x] =infinity Invalid argument 
| sin x * cosh y|_ >max_float Overflow 
| cos x * sinh y|_ >max_float Overflow 
See Also 


Appendix A, Critical Floating-Point Values 
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csqrt - Complex Square Root 


Interface 
F_COMPLEX csqrt (F_TYPE x, F_TYPE y) 


Description 


csqrt() computes the square root of a complex number, x + iy. The root is chosen 
so that the real part of csqrt(x,y) is greater than or equal to zero. 


Exceptions 
None. 
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cvt_ftof - Convert Between Supported Floating-Point Data Types 


Interface 


int cvt_ftof void *x, int x_type, void *y, int y_type, options 


Description 


Note 


This routine does not apply to OpenVMS Alpha. OpenVMS Alpha users 
should use the CVT$FTOF routine documented in the OpeVMS RTL 


Library (LIB$) Manual. 


cvt_ftof() converts a floating-point value from one data type to another. x points 
to the input value to be converted, and y points to the converted result. The 
conversion is subject to the options specified in the options (bit field) argument. 


x_type and y_type identify the data type of x and y as follows: 


Values for x_type and y_type 


Floating-Point Data Type 


CVT_VAX_F 

CVT_VAX_D 

CVT_VAX_G 

CVT_VAX_H 

CVT_IEEE_S 

CVT_IEEE_T 

CVT_IEEE_X 
CVT_BIG_ENDIAN_IEEE_S 
CVT_BIG_ENDIAN_IEEE_T 
CVT_BIG_ENDIAN_IEEE_X 
CVT_IBM_SHORT 
CVT_IBM_LONG 
CVT_CRAY_SINGLE 


VAX F Floating ( 4 bytes) 

VAX D Floating ( 8 bytes) 

VAX G Floating ( 8 bytes) 

VAX H Floating (16 bytes) 

IEEE Little Endian S Floating ( 4 bytes) 
IEEE Little Endian T Floating ( 8 bytes) 
IEEE Little Endian X Floating (16 bytes) 
IEEE Big Endian S Floating ( 4 bytes) 
IEEE Big Endian T Floating ( 8 bytes) 
IEEE Big Endian X Floating (16 bytes) 
IBM_Short_Floating ( 4 bytes) 
IBM_Long Floating ( 8 bytes) 

CRAY _Floating ( 8 bytes) 


Provide a zero (0) value to the options argument to select the default behavior 
or choose one or more options (status condition option, rounding options, 
"FORCE" options, CRAY and |1BM options) from the tables below as the options 
argument. Specify only the options that apply to your conversion. A conflicting 
or incompatible options argument will be reported as an error (CVT_INVALID_ 


OPTION). 
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Applicable 

Conversion Status Condition Option Description 

All CVT_REPORT_ALL Report all applicable status 

conditions as the default. The 
reporting of recoverable status 
conditions is disabled by default 
when this option is not used. 

Applicable 

Conversion Rounding Options Description 

All CVT_ROUND_TO_NEAREST The default rounding option for 

conversions to IEEE data types. 
This IEEE Std. 754 rounding 
mode results in the representable 
output value nearest to the 
infinitely precise result. If the two 
nearest representable values are 
equally near, the one with its least 
significant bit zero is the result. 

All CVT_BIASED_ROUNDING The default rounding option for 

conversions to non-IEEE data 

types. Performs "traditional" style 
rounding. This mode results in the 
reoresentable output value nearest 
to the infinitely precise result. If the 
two nearest representable values are 
equally near, the result is the value 
with the largest magnitude. 

All CVT_ROUND_TO_ZERO Round the output value toward zero 

(truncate). 

All CVT_ROUND_TO_POS Round the output value toward 

positive infinity. 

All CVT_ROUND_TO_NEG Round the output value toward 

negative infinity. 

Applicable 

Conversion "FORCE" Options Description 

All CVT_FORCE_ALL_SPECIAL_VALUES Apply all applicable "FORCE" 
options for the current conversion. 

IEEE CVT_FORCE_DENORM_TO ZERO! Force a denormalized |EEE output 
value to zero. 

IEEE CVT_FORCE_INF_TO MAX FLOAT? Force a positive IEEE infinite output 
value to +max_float and force a 
negative IEEE infinite output value 
to -max_float. 

IEEE or CVT_FORCE_INVALID_TO_ZERO2 Force an invalid IEEE NaN (not a 

VAX number) output value or a VAX ROP 


(reserved operand) output value to 
zero. 


1This option is valid only for conversions to IEEE output values. 
2This option is valid only for conversions to |EEE or VAX output values. 
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Returns 


Applicable 

Conversion Options for CRAY Format Conversion Description 

CRAY CVT_ALLOW_OVRFLW_RANGE_VALUES Allow an input/output 
exponent value > 60000 
(8). 

CRAY CVT_ALLOW_UDRFLW_RANGE_VALUES Allow an input/output 
exponent value < 20000 
(8). 

Applicable 

Conversion Option for IBM Format Conversion Description 

IBM CVT_ALLOW_UNNORMALIZED_VALUES Allow unnormalized input 


arguments. Allow an 
unnormalized output value 
for a small value that would 
normalize to zero. 


The return value is a bit field containing the condition codes raised by the 
function. cvt_ftof() returns CVT_NORMAL; otherwise, it sets one or more of the 
following recoverable and unrecoverable conditions. Use the following condition 
names to determine which conditions are set: 


Condition Name 


Condition (Always reported by default) 


CVT_INVALID_INPUT_TYPE 
CVT_INVALID_OUTPUT_TYPE 
CVT_INVALID_OPTION 


Invalid input type code. 
Invalid output type code. 


Invalid option argument. 


Condition Name 


Condition (Only reported if the CVT_ 
REPORT_ALL option is selected) 


CVT_RESULT_INFINITE 
CVT_RESULT_DENORMALIZED 
CVT_RESULT_OVERFLOW_RANGE 
CVT_RESULT_UNDERFLOW_RANGE 
CVT_RESULT_UNNORMALIZED 
CVT_RESULT_INVALID 


CVT_RESULT_OVERFLOW 
CVT_RESULT_UNDERFLOW 


Conversion produced an infinite result.? 
Conversion produced a denormalized result. 
Conversion yielded an exponent > 60000 (8).? 
Conversion yielded an exponent <20000 (8).2 
Conversion produced an unnormalized result.3 


Conversion result is either ROP (reserved operand), 
NaN (not a number), or closest equivalent. CRAY and 
IBM data types return 0.4 


Conversion resulted in overflow.* 


Conversion resulted in underflow.* 


1F or IEEE data type conversions. 
2F or CRAY data type conversions. 
3F or IBM data type conversions. 
4F or all data type conversions. 
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Condition (Only reported if the CVT_ 
Condition Name REPORT_ALL option is selected) 


CVT_RESULT_INEXACT Conversion resulted in a loss of precision.4 


4F or all data type conversions. 
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ANSI/IEEE Std 754-1985, IEEE Standard for Binary Floating-Point 
Arithmetic 


status = cvt_ftof( &big_x, CVT_BIG_ENDIAN_IEEE_T, 
&little_x, CVT_IEEE_T, 0 ); 


This example converts the value pointed to by big x, which is of type IEEE Big 
Endian T Floating, to the IEEE Little Endian T Floating data type. It stores the 
result in the location pointed to by little x. No conversion options are specified. 


status = cvt_ftof(&x, CVI_VAX_D, &y, CVI_IEEE_T, 
(CVT_FORCE_ALL_SPECIAL VALUES | CVI_REPORT_ALL) ); 


This example converts the value pointed to by x, which is of type VAX D Floating, 
to the IEEE Little Endian T Floating data type. It stores the result in the 
location pointed to by y. Any special IEEE values that would normally be 
generated will be removed. That is, NaN and Denormalized results will be 
returned as zero and infinite results will go to + max_float. In addition, all 
recordable status conditions will be reported. 
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drem - Remainder 


Interface 
F_TYPE drem (F_TYPE x, F_TYPE y) 


Description 


drem() returns the remainder r = x-n*y, where n = rint(x/y). Additionally, if 

| n-x/y| =1/2, then n is even. The remainder is computed exactly, and | r| is 
less than or equal to | y| /2. The drem() and remainder() functions are aliases 
of each other. 


Exceptions 


Exceptional Argument Routine Behavior 


x =infinity Invalid argument 


Note that rem(x,0) has value 0 and is not an exceptional case. 
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erf - Error Functions 


Interface 
F_TYPE erf (F_TYPE x) 
F_TYPE erfc (F_TYPE x) 


Description 


erf() returns the value of the error function. The definition of the erf() function is 
(2/sqrt(pi)) times the area under the curve exp(-t * t) between O and x. 


erfc() returns (1.0-erf(x)). 


Exceptions 
The erfc() function can result in an underflow as x gets large. 
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exp - Exponential 


Interface 
F_TYPE exp (F_TYPE x) 
F_TYPE expm1 (F_TYPE x) 
Description 
exp() computes the value of the exponential function, defined as e**x, where eis 
the constant used as a base for natural logarithms. 
expm1() computes exp(x)-1 accurately, even for tiny x. 
Exceptions 
Exceptional Argument Routine Behavior 
xX >In(max_float) Overflow 
x <In(min_float) Underflow 
See Also 


Appendix A, Critical Floating-Point Values 
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fabs - Absolute Value 


fabs - Absolute Value 


Interface 
F_TYPE fabs (F_TYPE x) 


Description 


fabs() computes the absolute value of x. 


Exceptions 


None. 
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finite - Check for Finite Value 


finite - Check for Finite Value 


Interface 
int finite (F_TYPE x) 


Description 
finite() returns the integer value 1 (true) or O (false). 
finite(x) = 1 when -infinity <x <-+infinity. 


finite(x) =O when | x| =infinity or x isa NaN. 


Exceptions 
None. 
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floor - Floor 


floor - Floor 


Interface 
F_TYPE floor (F_TYPE x) 


Description 


floor() returns the largest floating-point number of integral value less than or 
equal to x. 


Exceptions 


None. 
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fmod - Modulo Remainder 


Interface 
F_TYPE fmod (F_TYPE x, F_TYPE y) 


Description 


fmod() computes the floating-point remainder of x modulo y. It returns the 
remainder r =x-n*y, where n =trunc(x/y). The remainder is computed exactly. 


The result has the same sign as x and a magnitude less than the magnitude of y. 


Exceptions 


Exceptional Argument Routine Behavior 


x =infinity Invalid argument 


Note that fmod(x,0) has value O and is not an exceptional case. 
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fp_class - Classifies IEEE Floating-Point Values 


fp_class - Classifies IEEE Floating-Point Values 


Interface 


Description 


Exceptions 


See Also 


int fo_class (F_TYPE x) 


These routines determine the class of IEEE floating-point values. They return 
one of the constants in the file <fp_class.h>and never cause an exception, even 
for signaling NaNs. These routines implement the recommended function class(x) 
in the appendix of the IEEE Std 754. The constants in <fp class.h>refer to the 


following classes of values: 


Constant Class 

FP_SNAN Signaling NaN (Not-a-Number) 
FP_QNAN Quiet NaN (Not-a-Number) 
FP_POS_INF +H nfinity 

FP_NEG INF -I nfinity 

FP_POS NORM Positive normalized 


FP_NEG NORM 
FP_POS DENORM 
FP_NEG DENORM 


Negative normalized 
Positive denormalized 
Negative denormalized 


FP_POS ZERO +0.0 (positive zero) 
FP_NEG ZERO -0.0 (negative zero) 
None. 


ANSI/IEEE Std 754-1985, |EEE Standard for Binary Floating-Point Arithmetic 
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frexp - Convert to Fraction and Integral Power of 2 


Interface 
F_TYPE frexp (F_TYPE x, int *n) 


Description 


frexp() breaks a floating-point number into a normalized fraction and an integral 
power of 2. It stores the integer in the int object pointed to by the n parameter 
and returns the fraction part. 


Exceptions 
None. 
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hypot - Euclidean Distance 


Interface 
F_TYPE hypot (F_TYPE x, F_TYPE y) 

Description 
hypot() computes the length of the hypotenuse of a right triangle, where x and y 
represent the perpendicular sides of the triangle. 
hypot(x,y) is defined as the square root of (x**2 + y**2) and returns the same 
value as cabs(x,y). 

Exceptions 
Exceptional Argument Routine Behavior 
sqrt(x**2 + y**2) > max_float Overflow 

See Also 
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ilogb - Computes an Unbiased Exponent 


Interface 
int logo (F_TYPE x) 


Description 


ilogb(x) returns the unbiased exponent of x as an integer, (as if x were normalized 
>= 1.0 and <2.0) except: 


ilogb(NaN) is INT_MIN 
ilogb(inf) is INT_MAX 
logb(0) is INT_MIN 


There are no errors. The sign of x is ignored. 


Exceptions 


None. 
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isnan - Check for NaN Value 


Interface 
intisnan (F_TYPE x) 


Description 


isnan() returns 1 (true) if x is NaN (the IEEE floating-point reserved N ot-a- 
Number value) and 0 (false) otherwise. 


Exceptions 


None. 
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Idexp - Multiply by an Integral Power of 2 


Interface 
F_TYPE lIdexp (F_TYPE x, int n) 
Description 
Idexp() multiplies a floating-point number, x, by 2**n. 
Exceptions 
Exceptional Argument Routine Behavior 
| x*(2**n)| >max_float Overflow 
| x*(2**n) | <min_float Underflow 
See Also 
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Igamma - Computes the Logarithm of the gamma Function 


Interface 
F_TYPE lgamma (F_TYPE x) 


Description 


Igamma() returns the logarithm of the absolute value of gamma of x, or 

In(| G(x)| ), where G is the gamma function. The sign of gamma of x is returned 
in the external integer variable signgam as +1 or -1. The x parameter cannot be 
0 or a negative integer. 


gamma() returns the natural log of the gamma function and so is functionally 
equivalent to |gamma(). Because of this, gamma() is marked TO BE 
WITHDRAWN in the X/ Open Portability Guide, Revision 4 (XPG4). 


Exceptions 


Exceptional Argument Routine Behavior 


| x| =infinity Invalid argument 
x =0, -1, -2, -3, ... Invalid argument 
| x| >lgamma_max_float Overflow 


See Also 
Appendix A, Critical Floating-Point Values 
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log - Logarithm Functions 


log - Logarithm Functions 


Interface 
F_TYPE In (F_TYPE x) 
F_TYPE log2 (F_TYPE x) 
F_TYPE logi0 (F_TYPE x) 
F_TYPE logip (F_TYPE y) 


Description 
In() computes the natural (base e) logarithm of x. 
log2() computes the base 2 logarithm of x. 
log10() computes the common (base 10) logarithm of x. 
loglp( ) computes In(1+y) accurately, even for tiny y. 


Exceptions 


Exceptional Argument Routine Behavior 
x <0 Invalid argument 
x =0 Overflow 

1+y <0 Invalid argument 
1+y =0 Overflow 


CPML-46 CPML Routines 


CPML Routines 
logb - Radix-independent Exponent 


logb - Radix-independent Exponent 


Interface 
F_TYPE logb (F_TYPE x) 


Description 


logb() returns a signed integer converted to double-precision floating-point and so 
chosen that 1 <= | x| /2**n <2 unless x =Oor | x| =infinity. 


IEEE Std 754 defines logb(+infinity) = +infinity and logb(0) =-infinity. The latter 
is required to signal division by zero. 


Exceptions 


Exceptional Argument Routine Behavior 


x =0 Invalid argument 
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modf - Return the Fractional Part and Integer Part of a Floating-Point 
Number 


Interface 
F_TYPE modf (F_TYPE x, F_TYPE *n) 


Description 


modf() splits a floating-point number x into a fractional part f and an integer 
part i such that | f| <1.0 and (f +i) =x. Both f andi have the same sign as x. 
modf() returns f and stores i into the location pointed to by n. 


Exceptions 


None. 
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nextafter - Next Machine Number After 


Interface 
F_TYPE nextafter (F_TYPE x, F_TYPE y) 
Description 
nextafter() returns the machine-representable number next to x in the direction 
y. 
Exceptions 
Exceptional Argument Routine Behavior 
X =max_float and y = +infinity Overflow 
X =-max_float and y =-infinity Overflow 
xX =min_float and y is less than or equal to 0 Underflow 
xX =-min_float and y is greater than or equal to 0 Underflow 
See Also 


ANSI/IEEE Std 754-1985, [EEE Standard for Binary Floating-Point Arithmetic 
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nint - Round to the Nearest Integer 


nint - Round to the Nearest Integer 


Interface 
F_TYPE nint (F_TYPE x) 


Description 


nint() returns the nearest integral value to x, except halfway cases are rounded to 
the integral value larger in magnitude. This function corresponds to the Fortran 
generic intrinsic function nint(). 


Exceptions 
None. 
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pow - Raise the Base to a Floating-Point Exponent 


Interface 
F_TYPE pow (F_TYPE x, F_TYPE y) 
Description 
pow() raises a floating-point base x to a floating-point exponent y. The value of 
pow(x,y) is computed as e**(y In(x)) for positive x. If x is O or negative, see your 
language reference manual. 
Passing a NaN input value to pow() produces a NaN result for nonzero values of 
y. For pow(NaN,0), see your language reference manual. 
Exceptions 
Exceptional Argument Routine Behavior 
y In(x) >In(max_float) Overflow 
y In(x) <In(min_float) Underflow 
Fortran-Exceptional Argument Routine Behavior 
x <0 Invalid argument 
x =O and y <0 Invalid argument 
x =O and y =0 Invalid argument 
x =4infinity and y =0 Invalid argument 
xX =land | y| =infinity Invalid argument 
ANSI C-Exceptional Argument Routine Behavior 
| x| =land | y| =infinity Invalid argument 
x <0 and y is not integral Invalid argument 
See Also 
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random - Random Number Generator, Uniformly Distributed 


Interface 
F_TYPE random (int *n) 


Description 


random() is a general random number generator. The argument to the random 
function is an integer passed by reference. There are no restrictions on the input 
argument, although it should be initialized to different values on separate runs in 
order to obtain different random sequences. This function must be called again to 
obtain the next pseudo random number. The argument is updated automaticall y. 


The result is a floating-point number that is uniformly distributed in the interval 
(0.0,1.0). 


Exceptions 
None. 
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remainder - Remainder 


Interface 
F_TYPE remainder (F_TYPE x, F_TYPE y) 


Description 


remainder() returns the remainder r =x-n*y, where n =rint(x/y). Additionally, if 
| n-x/y| =1/2, then n is even. Consequently, the remainder is computed exactly, 
and | r| is less than or equal to | y| /2. The drem() and remainder() functions 
are aliases of each other. 


Exceptions 


Exceptional Argument Routine Behavior 


x =infinity Invalid argument 


Note that rem(x,0) has value 0 and is not an exceptional case. 
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rint - Return the Nearest Integral Value 


Interface 
F_TYPE rint (F_TYPE x) 


Description 


rint() rounds x to an integral value according to the current IEEE rounding 
direction specified by the user. 


Exceptions 


None. 
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scalb - Exponent Adjustment 


Interface 
F_TYPE scalb (F_TYPE x, F_TYPE y) 
Description 
scalb() = x*(2**y) computed, for integer-valued floating point number y. 
Exceptions 
Exceptional Argument Routine Behavior 
x*(2**y) > max_float Overflow 
x*(2**y) <min_float Underflow 
x=0, y=infinity Invalid argument 
x=infinity, y=infinity Invalid argument 
See Also 
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sin - Sine of Angle 


sin - Sine of Angle 


Interface 


Description 


Exceptions 


See Also 


F_TYPE sin (F_TYPE x) 
F_TYPE sind (F_TYPE x) 


sin() computes the sine of x, measured in radians. 


sind() computes the sine of x, measured in degrees. 


Exceptional Argument Routine Behavior 
| x] =infinity Invalid argument 
(sind) | x] <(180/pi) * min_float Underflow 
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sincos - Sine and Cosine of Angle 


sincos - Sine and Cosine of Angle 


Interface 


Description 


Exceptions 


F_COMPLEX sincos (F_TYPE x) 
F_COMPLExX sincosd (F_TYPE x) 


sincos() computes both the sine and cosine of x, measured in radians. 
sincosd() computes both the sine and cosine of x, measured in degrees. 
sincos(x) is defined as (sin x +iCos y). 


Exceptional Argument Routine Behavior 
| x| =infinity Invalid argument 
(sind) | x] <(180/pi) * min_float Underflow 
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sinh - Hyperbolic Sine 


Interface 
F_TYPE sinh (F_TYPE x) 
Description 
sinh() computes the hyperbolic sine of x. 
sinh(x) is defined as (exp(x)-exp(-x))/2. 
Exceptions 
Exceptional Argument Routine Behavior 
| x| >In(2 * max_float) Overflow 
See Also 
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sinhcosh - Hyperbolic Sine and Cosine 


sinhcosh - Hyperbolic Sine and Cosine 


Interface 
F_COMPLEX sinhcosh (F_TYPE x) 


Description 
sinhcosh( ) computes both the hyperbolic sine and hyperbolic cosine of x. 


sinhcosh(x) is defined as (sinh x +icosh x). 


Exceptions 


Exceptional Argument Routine Behavior 


| x] >In(2 * max_float) Overflow 


See Also 
Appendix A, Critical Floating-Point Values 
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sqrt - Square Root 


sqrt - Square Root 


Interface 
F_TYPE sqrt (F_TYPE x) 


Description 
sqrt() computes the rounded square root of x. 


For platforms supporting a signed zero, sqrt(-0) =0. 


Exceptions 


Exceptional Argument Routine Behavior 


x <0 Invalid argument 
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tan - Tangent of Angle 


Interface 
F_TYPE tan (F_TYPE x) 
F_TYPE tand (F_TYPE x) 
Description 
tan() computes the tangent of x, measured in radians. 
tand() computes the tangent of x, measured in degrees. 
Exceptions 
Exceptional Argument Routine Behavior 
| x] =infinity Invalid argument 
(tand) | x] <(180/pi) * min_float Underflow 
(tand) x =(2n+1) * 90 Overflow 
See Also 
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tanh - Hyperbolic Tangent 


tanh - Hyperbolic Tangent 


Interface 
F_TYPE tanh (F_TYPE x) 


Description 
tanh() computes the hyperbolic tangent of x. 
tanh(x) is defined as (exp(x)-exp(-x))/(exp(x) + exp(-x)). 


Exceptions 


None. 
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trunc - Truncation 


Interface 
F_TYPE trunc (F_TYPE x) 


Description 


trunc() truncates x to an integral value. 


Exceptions 


None. 
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unordered - Check for x Unordered with Respect to y 


Interface 
int unordered (F_TYPE x, F_TYPE y) 


Description 


unordered(x,y) returns the value 1 (true) if x, y, or both area NaN and returns 
the value 0 (false) otherwise. 


Exceptions 


None. 


CPML-64 CPML Routines 


A 


Critical Floating-Point Values 


Table A-1 contains the hexadecimal and decimal boundary values used in CPML 
calculations and exception checking. 


Table A-—1 Hexadecimal and Decimal Boundary Values 


Data 

Type Value for: max_float 

F Hexadecimal: FF FF 7FFF 

G Hexadecimal: FF FFFFFFFFFF 7FFF 

S Hexadecimal: 7F 7F FF FF 

T Hexadecimal: 7FEFFFFFFFFFFFFF 

x Hexadecimal: 7FFEFFFFFFFFFFFFFFFFFFFFFFFFFFFF 
F Decimal: 1.701411e38 

G Decimal: 8.988465674311579e307 

S Decimal: 3.402823e38 

T Decimal: 1.797693134862316e308 

x Decimal: 1.189731495357231765085759326628007016196477e4932 


Data 

Type Value for: min_float 

F Hexadecimal: 00000080 

G Hexadecimal: 0000000000000010 

S Hexadecimal: 00000001 

T Hexadecimal: 0000000000000001 

x Hexadecimal: 00000000000000000000000000000001 
F Decimal: 2.9387359e-39 

G Decimal: 5.562684646268003e-309 

S) Decimal: 1.4012985e-45 

T Decimal: 4.940656458412465e-324 

x Decimal: 6.4751751194380251109244389582276465524996e-4966 


(continued on next page) 
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Table A—1 (Cont.) Hexadecimal and Decimal Boundary Values 


Data 

Type Value for: In(max_float) 

F Hexadecimal: OF 3443B0 

G Hexadecimal: 7B616E 3A28B 740A6 

S Hexadecimal: 42B17218 

T Hexadecimal: 40862E 42F EFA39E F 

x Hexadecimal: 400C62E 42F EFA39E F 35793C7673007E 6 


Decimal: 88.029692 

Decimal: 709.0895657128241 

Decimal: 88.7228391 

Decimal: 709.7827128933840 

Decimal: 11356.5234062941439494919310779707648912527 


<x AWNQa Tm 


Type Value for: In(min_float) 
F Hexadecimal: 7218C3B1 


G Hexadecimal: 39EFFEFA2E 42COA6 
S Hexadecimal: C2CE8EDO 
T 
X 


Hexadecimal: C0874385446D71C3 
Hexadecimal: C00C6546282207802C 89D 24D65E 96274 


Decimal: -88.72284 
Decimal: -709.7827128933840 
Decimal: -103.2789 
Decimal: -744.4400719213813 
Decimal: -11432.7695961557379335278266113311643138373 
(continued on next page) 
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Table A-1 (Cont.) Hexadecimal and Decimal Boundary Values 


Data 

Type Value for: In(2 * max_float) 

F Hexadecimal: 721843B1 

G Hexadecimal: 39E FFEFA2E 4240A6 

S Hexadecimal: 42B2D4F C 

T Hexadecimal: 408633CE 8F B9F 87E 

X Hexadecimal: 400C62E 9BB 80635D81D36125B 64DA4A6 
F Decimal: 88.72284 

G Decimal: 709.7827128933840 

S) Decimal: 89.41599 

T Decimal: 710.4758600739439 

x Decimal: 11357.2165534747038948013483100922230678208 
Data 

Type Value for: (180/pi) * min_float 

F Hexadecimal: 2EE 10365 

G Hexadecimal: C1F 81A63A5DCO006C 

S Hexadecimal: 00000039 

T Hexadecimal: 0000000000000039 

x Hexadecimal: 00000000000000000000000000000039 

F Decimal: 1.683772e-37 

G Decimal: 3.187183529933798e-307 

S Decimal: 8.028849e-44 

T Decimal: 2.830787630910868e-322 

x Decimal: 3.71000205951917569316937757202433432154392e-4964 


(continued on next page) 
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Table A-—1 (Cont.) Hexadecimal and Decimal Boundary Values 


Data 

Type Value for: Igamma_max_float 

F Hexadecimal: 50F 97CC6 

G Hexadecimal: F55F C5015ABD7F 67 

S Hexadecimal: 7BC650F 9 

T Hexadecimal: 7F 4754BDC501F 55F 

X Hexadecimal: 7FF171AA9917F FF BD7EA44AE 6D203DF 6 
F Decimal: 2.0594342e36 

G Decimal: 1.2812545499066958e305 

S Decimal: 2.0594342e36 

T Decimal: 1.2812545499066958e305 

X Decimal: 1.0485738685148938358098967157129705040168e4928 


A-4 Critical Floating-Point Values 


B 


CPML Entry-Point Names 


Each entry-point name in Table B-1 is unique and corresponds to data-type 
specific calculations in a CPML routine. For example, the acos function has 
five entry-point-names for the OpenVMS Alpha operating system. Because 
five floating-point data types are available, five acos routines are provided: 
math$acos_s, math$acos t, math$acos f, math$acos_g, and math$acos_ x. Use 
the entry-point name that corresponds to your input argument data type. 


Table B—1_ Entry-Point Names for CPML Platforms 


Entry-Point Names 


Generic 

Function Data Type Compaq Tru64 UNIX 

Name Required OpenVMS Alpha Alpha 

acos S_ FLOAT math$acos_s acosf 
T_FLOAT math$acos_t acos 
X_FLOAT math$acos_x acosl 
F_FLOAT math$acos_f 
G_FLOAT math$acos_g 

acosd S_ FLOAT math$acosd_s acosdf 
T_FLOAT math$acosd_t acosd 
X_FLOAT math$acosd_x acosdl 
F_FLOAT math$acosd_f 
G_ FLOAT math$acosd_g 

acosh S_ FLOAT math$acosh_s acoshf 
T_FLOAT math$acosh_t acosh 
X_FLOAT math$acosh_x acoshl 
F_FLOAT math$acosh_f 
G_ FLOAT math$acosh_g 


(continued on next page) 
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Table B—1 (Cont.) Entry-Point Names for CPML Platforms 


Entry-Point Names 


Generic 

Function Data Type Compaq Tru64 UNIX 

Name Required OpenVMS Alpha Alpha 

asin S_ FLOAT math$asin_s asinf 
T_FLOAT math$asin_t asin 
X_FLOAT math $asin_x asin 
F_FLOAT math$asin_f 
G_ FLOAT math$asin_g 

asind S_ FLOAT math$asind_s asindf 
T_FLOAT math$asind_t asind 
X_FLOAT math$asind_x asindl 
F_FLOAT math$asind_f 
G_ FLOAT math$asind_g 

asinh S_FLOAT math$asinh_s asinhf 
T_FLOAT math$asinh_t asinh 
X_FLOAT math$asinh_x asinhl 
F_FLOAT math$asinh_f 
G_ FLOAT math$asinh_g 

atan S_ FLOAT math$atan_s atanf 
T_FLOAT math$atan_t atan 
X_FLOAT math$atan_x atanl 
F_FLOAT math$atan_f 
G_FLOAT math$atan_g 

atan2 S_ FLOAT math$atan2_s atan2f 
T_FLOAT math$atan2 t atan2 
X_FLOAT math$atan2_x atan2l 
F_FLOAT math$atan2_f 
G_FLOAT math$atan2_g 

atand S_ FLOAT math $atand_s atandf 
T_FLOAT math$atand_t atand 
X_FLOAT math$atand_x atandl 
F_FLOAT math$atand_f 
G_FLOAT math$atand_g 
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Table B—1 (Cont.) Entry-Point Names for CPML Platforms 


Entry-Point Names 


Generic 

Function Data Type Compaq Tru64 UNIX 

Name Required OpenVMS Alpha Alpha 

atand2 S_ FLOAT math$atand2_s atand2f 
T_FLOAT math$atand2_t atand2 
X_FLOAT math$atand2_x atand2l 
F_FLOAT math$atand2_f 
G_ FLOAT math$atand2_g 

atanh S_ FLOAT math$atanh_s atanhf 
T_FLOAT math$atanh_t atanh 
X_FLOAT math$atanh_x atanhl 
F_FLOAT math$atanh_f 
G_ FLOAT math$atanh_g 

cabs S_ FLOAT math$hypot_s cabsf 
T_FLOAT math$hypot_t cabs 
X_FLOAT math$hypot_x cabsl 
F_FLOAT math$hypot_f 
G_ FLOAT math$hypot_g 

cort S_ FLOAT math$cbrt_s cortf 
T_FLOAT math$cbrt_t cort 
X_FLOAT math$cbrt_x cbrtl 
F_FLOAT math$cbrt_f 
G_ FLOAT math$cbrt_g 

ccos S_ FLOAT math$ccos_s ccosf 
T_FLOAT math$ccos_t ccos 
X_FLOAT math$ccos_x ccosl 
F_FLOAT math$ccos_f 
G_ FLOAT math$ccos_g 
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Table B—1 (Cont.) Entry-Point Names for CPML Platforms 


Entry-Point Names 


Generic 

Function Data Type Compaq Tru64 UNIX 

Name Required OpenVMS Alpha Alpha 

cdiv S_ FLOAT math$cdiv_s cdi vf 
T_FLOAT math$cdiv_t cdiv 
X_FLOAT math $cdiv_x cdivl 
F_FLOAT math$cdiv_f 
G_FLOAT math$cdiv_g 

ceil S_ FLOAT math $ceil_s ceilf 
T_FLOAT math$ceil_t ceil 
X_FLOAT math$ceil_x ceill 
F_FLOAT math$ceil_f 
G_FLOAT math$ceil_g 

cexp S_FLOAT math$cexp_s cexpf 
T_FLOAT math$cexp_t cexp 
X_FLOAT math$cexp_x cexpl 
F_FLOAT math$cexp_f 
G_FLOAT math$cexp_g 

dog S_ FLOAT math$clog_s clogf 
T_FLOAT math$clog_t clog 
X_FLOAT math$clog_x clogl 
F_FLOAT math$clog_f 
G_FLOAT math$clog_g 

cmul S_ FLOAT math$cmul_s cmulf 
T_FLOAT math$cmul_t cmul 
X_FLOAT math$cmul_x cmull 
F_FLOAT math$cmul_f 
G_FLOAT math$cmul_g 

copysign S_ FLOAT math$copysign_s copysi gnf 
T_FLOAT math$copysign_t copysign 
X_FLOAT math$copysign_x copysign 
F_FLOAT math$copysign_f 
G_FLOAT math$copysign_g 
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Entry-Point Names 


Generic 

Function Data Type Compaq Tru64 UNIX 

Name Required OpenVMS Alpha Alpha 

cos S_ FLOAT math$cos_s cosf 
T_FLOAT math$cos_t cos 
X_FLOAT math$cos_x cos| 
F_FLOAT math$cos_f 
G_ FLOAT math$cos_g 

cosd S_ FLOAT math$cosd_s cosdf 
T_FLOAT math$cosd_t cosd 
X_FLOAT math$cosd_x cosdl 
F_FLOAT math$cosd_f 
G_ FLOAT math$cosd_g 

cosh S_ FLOAT math$cosh_s coshf 
T_FLOAT math$cosh_t cosh 
X_FLOAT math$cosh_x cosh 
F_FLOAT math$cosh_f 
G_ FLOAT math$cosh_g 

cot S_ FLOAT math$cot_s cotf 
T_FLOAT math$cot_t cot 
X_FLOAT math$cot_x cotl 
F_FLOAT math$cot_f 
G_ FLOAT math$cot_g 

cotd S_ FLOAT math$cotd_s cotdf 
T_FLOAT math$cotd_t cotd 
X_FLOAT math$cotd_x cotdl 
F_FLOAT math$cotd_f 
G_ FLOAT math$cotd_g 
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Entry-Point Names 


Generic 

Function Data Type Compaq Tru64 UNIX 

Name Required OpenVMS Alpha Alpha 

cpow S_ FLOAT math$cpow_s cpowf 
T_FLOAT math$cpow_t cpow 
X_FLOAT math$cpow_x cpowl 
F_FLOAT math$cpow_f 
G_FLOAT math$cpow_g 

csin S_ FLOAT math $csin_s csinf 
T_FLOAT math$csin_t csin 
X_FLOAT math $csin_x csinl 
F_FLOAT math$csin_f 
G_ FLOAT math$csin_g 

csqrt S_ FLOAT math$csqrt_s csqrtf 
T_FLOAT math$csqrt_t csqrt 
X_FLOAT math$csqrt_x csqrtl 
F_FLOAT math$csqrt_f 
G_ FLOAT math$csqrt_g 

cvt_ftof All supported cvt_ftof 
types 

drem S_ FLOAT math$rem_s dremf 
T_FLOAT math$rem_t drem 
X_FLOAT math$rem_x dreml 
F_FLOAT math$rem_f 
G_FLOAT math$rem_g 

erf S_ FLOAT math$erf_s erff 
T_FLOAT math$erf_t erf 
X_FLOAT math $erf_x erfl 
F_FLOAT math $erf_f 
G_FLOAT math $erf_g 
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Entry-Point Names 


Generic 

Function Data Type Compaq Tru64 UNIX 

Name Required OpenVMS Alpha Alpha 

erfc S_ FLOAT math$erfc_s erfcf 
T_FLOAT math$erfc_t erfc 
X_FLOAT math$erfc_x erfdl 
F_FLOAT math$erfc_f 
G_ FLOAT math$erfc_g 

exp S_ FLOAT math$exp_s expf 
T_FLOAT math$exp_t exp 
X_FLOAT math$exp_x expl 
F_FLOAT math$exp_f 
G_ FLOAT math$exp_g 

expm1 S_ FLOAT math$expm1_s expmlf 
T_FLOAT math$expm1_t expm1 
X_FLOAT math$expm1_x expm1l 
F_FLOAT math$expm1_f 
G_ FLOAT math$expm1_g 

fabs S_ FLOAT math$fabs_s fabsf 
T_FLOAT math$fabs t fabs 
X_FLOAT math$fabs_x fabsl 
F_FLOAT math$fabs_f 
G_ FLOAT math$fabs_g 

finite S_ FLOAT math$finite_s finitef 
T_FLOAT math$finite_t finite 
X_FLOAT math$finite_x finitel 
F_FLOAT math$finite_f 
G_ FLOAT math$finite_g 
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Entry-Point Names 


Generic 

Function Data Type Compaq Tru64 UNIX 

Name Required OpenVMS Alpha Alpha 

floor S_ FLOAT math $floor_s floor 
T_FLOAT math$floor_t floor 
X_FLOAT math $floor_x floor| 
F_FLOAT math$floor_f 
G_FLOAT math$floor_g 

fmod S_ FLOAT math$mod_s fmodf 
T_FLOAT math$mod_t fmod 
X_FLOAT math$mod_x fmodl 
F_FLOAT math $mod_f 
G_FLOAT math$mod_g 

fp_class S_FLOAT math$fp_class_s fp_classf 
T_FLOAT math$fp_class t fp_class 
X_FLOAT math$fp_class_x fp_classl 
F_FLOAT math$fp_class_f 
G_FLOAT math$fp_class_g 

frexp S_ FLOAT math$frexp_s frexpf 
T_FLOAT math $frexp_t frexp 
X_FLOAT math $frexp_x frexpl 
F_FLOAT math $frexp_f 
G_FLOAT math$frexp_g 

hypot S_ FLOAT math$hypot_s hypotf 
T_FLOAT math$hypot_t hypot 
X_FLOAT math$hypot_x hypotl 
F_FLOAT math$hypot_f 
G_FLOAT math$hypot_g 

ilogb S_ FLOAT math$ilogb_s ilogbf 
T_FLOAT math$ilogb_t ilogb 
X_FLOAT math$ilogb_x ilogbl 
F_FLOAT math$ilogb_f 
G_ FLOAT math$ilogb_g 
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Entry-Point Names 


Generic 

Function Data Type Compaq Tru64 UNIX 

Name Required OpenVMS Alpha Alpha 

isnan S_ FLOAT math$isnan_s isnanf 
T_FLOAT math$isnan_t isnan 
X_FLOAT math$isnan_x isnanl 
F_FLOAT math$isnan_f 
G_ FLOAT math$isnan_g 

jo S_FLOAT math$j0_s jOf 
T_FLOAT math$j0_t jo 
X_FLOAT math$jO_x jol 
F_FLOAT math$jO_f 
G_ FLOAT math$j0O_g 

jl S_FLOAT math$j1_s jif 
T_FLOAT math$j1_t jl 
X_FLOAT math$j1_x jul 
F_FLOAT math$j1_f 
G_ FLOAT math$j1_g 

jn S_ FLOAT math$jn_s jnf 
T_FLOAT math$jn_t jn 
X_FLOAT math$jn_x jnl 
F_FLOAT math$jn_f 
G_ FLOAT math$jn_g 

Idexp S_ FLOAT math$ldexp_s Idexpf 
T_FLOAT math$ldexp_t Idexp 
X_FLOAT math$ldexp_x Idexpl 
F_FLOAT math$ldexp_f 
G_ FLOAT math$ldexp_g 
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Generic 

Function Data Type Compaq Tru64 UNIX 

Name Required OpenVMS Alpha Alpha 

Igamma S_ FLOAT math$lgamma_s Igammaf 
T_FLOAT math$lgamma_t Igamma 
X_FLOAT math$lgamma_x Igammal 
F_FLOAT math$lgamma_f 
G_FLOAT math$lgamma_g 

In S_ FLOAT math$In_s loof 
T_FLOAT math$In_t log 
X_FLOAT math$In_x logl 
F_FLOAT math$lIn_f 
G_FLOAT math$In_g 

log2 S_FLOAT math$log2_s log2f 
T_FLOAT math$log2_t log2 
X_FLOAT math$log2_x log2| 
F_FLOAT math$log2_f 
G_FLOAT math$log2_g 

log10 S_ FLOAT math$log10_s loglOf 
T_FLOAT math$log10_t log10 
X_FLOAT math$log10_x log10l 
F_FLOAT math$log10_f 
G_ FLOAT math$logl10_g 

loglp S_ FLOAT math$loglp_s loglpf 
T_FLOAT math$loglp_t loglp 
X_FLOAT math$loglp_x loglpl 
F_FLOAT math$loglp_f 
G_FLOAT math$loglp_g 

logb S_ FLOAT math$logb_s logbf 
T_FLOAT math$logb_t logb 
X_FLOAT math$logb_x logbl 
F_FLOAT math$logb_f 
G_FLOAT math$logb_g 


B-10 CPML Entry-Point Names 


(continued on next page) 


CPML Entry-Point Names 


Table B—1 (Cont.) Entry-Point Names for CPML Platforms 


Entry-Point Names 


Generic 

Function Data Type Compaq Tru64 UNIX 

Name Required OpenVMS Alpha Alpha 

modf S_ FLOAT math$modf_s modff 
T_FLOAT math$modf_t modf 
X_FLOAT math$modf_x modfl 
F_FLOAT math$modf_f 
G_ FLOAT math$modf_g 

nextafter S_ FLOAT math$nextafter_s nextafterf 
T_FLOAT math$nextafter_t nextafter 
X_FLOAT math$nextafter_x nextafter| 
F_FLOAT math$nextafter_f 
G_ FLOAT math$nextafter_g 

nint S_ FLOAT math$nint_s nintf 
T_FLOAT math$nint_t nint 
X_FLOAT math$nint_x nintl 
F_FLOAT math$nint_f 
G_ FLOAT math$nint_g 

pow S_ FLOAT math$pow_ss powf 
T_FLOAT math$pow_tt pow 
X_FLOAT math$pow_xx powl 
F_FLOAT math$pow_ff 
G_ FLOAT math$pow_gg 

random S_ FLOAT math$random_l_s 
T_FLOAT 
X_FLOAT 
F_FLOAT math$random_|_f 
G FLOAT 
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Generic 

Function Data Type Compaq Tru64 UNIX 

Name Required OpenVMS Alpha Alpha 

remainder S_ FLOAT math$rem_s remainderf 
T_FLOAT math$rem_t remainder 
X_FLOAT math$rem_x remainder! 
F_FLOAT math$rem_f 
G_FLOAT math$rem_g 

rint S_ FLOAT math$rint_s rintf 
T_FLOAT math$rint_t rint 
X_FLOAT math$rint_x rintl 
F_FLOAT math$rint_f 
G_FLOAT math$rint_g 

scalb S_ FLOAT math$scalb_s scalbf 
T_FLOAT math$scalb_t scalb 
X_FLOAT math$scalb_x scalbl 
F_FLOAT math$scalb_f 
G_ FLOAT math$scalb_g 

sin S_ FLOAT math$sin_s sinf 
T_FLOAT math$sin_t sin 
X_FLOAT math$sin_x sinl 
F_FLOAT math$sin_f 
G_FLOAT math$sin_g 

sincos S_ FLOAT math $sincos_s sincosf 
T_FLOAT math $sincos_t sincos 
X_FLOAT math$sincos_x sincosl 
F_FLOAT math$sincos_f 
G_FLOAT math$sincos_g 

sincosd S_ FLOAT math $sincosd_s sincosdf 
T_FLOAT math$sincosd_t sincosd 
X_FLOAT math $sincosd_x sincosdl 
F_FLOAT math$sincosd_f 
G_ FLOAT math$sincosd_g 
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Generic 
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Name Required OpenVMS Alpha Alpha 

sind S_ FLOAT math$sind_s sindf 
T_FLOAT math$sind_t sind 
X_FLOAT math$sind_x sindl 
F_FLOAT math$sind_f 
G_ FLOAT math$sind_g 

sinh S_ FLOAT math$sinh_s sinhf 
T_FLOAT math$sinh_t sinh 
X_FLOAT math$sinh_x sinhl 
F_FLOAT math$sinh_f 
G_ FLOAT math$sinh_g 

sinhcosh S_ FLOAT math$sinhcosh_s sinhcoshf 
T_FLOAT math$sinhcosh_t sinhcosh 
X_FLOAT math$sinhcosh_x sinhcoshl 
F_FLOAT math$sinhcosh_f 
G_ FLOAT math$sinhcosh_g 

sqrt S_ FLOAT math$sqrt_s sqrtf 
T_FLOAT math$sqrt_t sqrt 
X_FLOAT math$sqrt_x sqrtl 
F_FLOAT math$saqrt_f 
G_ FLOAT math$sqrt_g 

tan S_ FLOAT math$tan_s tanf 
T_FLOAT math$tan_t tan 
X_FLOAT math$tan_x tanl 
F_FLOAT math$tan_f 
G_ FLOAT math$tan_g 


(continued on next page) 


CPML Entry-Point Names B-13 


CPML Entry-Point Names 


Table B—1 (Cont.) Entry-Point Names for CPML Platforms 


Entry-Point Names 


Generic 

Function Data Type Compaq Tru64 UNIX 

Name Required OpenVMS Alpha Alpha 

tand S_ FLOAT math$tand_s tandf 
T_FLOAT math$tand_t tand 
X_FLOAT math$tand_x tandl 
F_FLOAT math$tand_f 
G_FLOAT math$tand_g 

tanh S_ FLOAT math$tanh_s tanhf 
T_FLOAT math$tanh_t tanh 
X_FLOAT math$tanh_x tanhl 
F_FLOAT math$tanh_f 
G_FLOAT math$tanh_g 

trunc S_FLOAT math$trunc_s truncf 
T_FLOAT math$trunc_t trunc 
X_FLOAT math$trunc_x truncl 
F_FLOAT math$trunc_f 
G_FLOAT math$trunc_g 

unordered S_ FLOAT math $unordered_s unorderedf 
T_FLOAT math$unordered_t unordered 
X_FLOAT math $unordered_x unordered 
F_FLOAT math$unordered_f 
G_FLOAT math$unordered_g 

yO S_FLOAT math$y0_s yOf 
T_FLOAT math$y0 t yO 
X_FLOAT math$y0_x yol 
F_FLOAT math$y0_f 
G_FLOAT math$y0_g 

yl S_ FLOAT math$y1_s ylf 
T_FLOAT math$y1 t yl 
X_FLOAT math$y1_x yll 
F_FLOAT math$y1_f 
G_FLOAT math$yl1_g 
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Generic 
Function Data Type Compaq Tru64 UNIX 
Name Required OpenVMS Alpha Alpha 
yn S_ FLOAT math$yn_s ynf 
T_FLOAT math$yn_t yn 
X_FLOAT math$yn_x ynl 
F_FLOAT math$yn_f 
G_FLOAT math$yn_g 


CPML Entry-Point Names B-15 


Glossary 


This glossary defines mathematical terms and symbolic names used in this 
manual. 


complex number 
See F_ COMPLEX. 


denormalized number 
A floating-point number with a value very close to zero. 


domain error 

An exception condition resulting from passing an argument whose value is 
outside the range of permissible values. 

exceptional argument 

Any argument value passed to a CPML routine that does not return a meaningful 
result, or an argument defined differently for different platforms. 

F_COMPLEX 


A complex number identifier. F_COMPLEX indicates that a given routine returns 
two different values of the same floating-point data type. See Table 1-2 for more 
information. 


F_TYPE 


A floating-point number identifier. F_TYPE is used when it is not necessary to 
distinguish between the floating types. See Table 1-1 for more information. 


floating-point number 
See F_TYPE. 


HUGE_RESULT 

For VAX data types, HUGE_RESULT =max_float. 
For IEEE data types, HUGE RESULT =infinity. 
INV_RESULT 

For VAX data types, INV_RESULT =0. 

For IEEE data types, INV_RESULT =NaN. 


invalid argument 
See domain error. 


Glossary-1 


max_float 


The largest finite number representable in the floating-point data types. See 
Appendix A for more information on max_float values. 


min_float 


The smallest positive normalized nonzero number representable in the 
floating-point data types. See Appendix A for more information on min_float 
values. 


NaN 


A floating-point value that is said to be “not a number” and contains an 
indeterminate quantity. 


overflow 


An exception condition caused by passing a floating-point value that is larger 
than the highest valid floating-point value. See max_float for additional 
information. 


range error 


An exception condition that occurs when a mathematically valid argument 
results in a function value that exceeds the range of representable values for 
floating-point data types. 


underflow 


An exception condition caused by passing a floating-point value that is lower than 
the lowest valid floating-point value. See min_float for additional information. 


Glossary-2 


A 


32-bit IEEE single-precision complex number, 1-3 
32-bit IEEE single-precision number, 1-2 
Absolute value, CPML-35 
32-bit VAX single-precision complex number, 1-3 
acosdf routine, B-1 
acosd routine, CPML-5, B-1 
acosf routine, B-1 
acoshf routine, B-1 
acosh routine, CPML-6, B-1 
acos routine, CPML-5, B-1 
32-bit VAX single-precision number, 1-2 
64-bit IEEE double-precision complex number, 
1-3 
64-bit IEEE double-precision number, 1-2 
64-bit VAX double-precision complex number, 1-3 
64-bit VAX double-precision number, 1-2 
128-bit IEEE extended-precision complex number, 
1-3 
128-bit IEEE extended-precision number, 1-2 
Arc cosine of angle, CPML-5 
hyperbolic, CPML-6 
Arc sine of angley CPML-7 
Arc tangent of angle 
hyperbolic, CPML-11 
with one argument, CPML-9 
with two arguments, CPML-10 
asindf routine, B-2 
asind routine, CPML-7, B-2 
asinf routine, B-2 
asinhf routine, B-2 
asinh routine, CPML-8, B-2 
asin routine, CPML-7, B-2 
atan2f routine, B-2 
atan2 routine, 2-3, CPML-10, B-2 
atand2f routine, B-3 
atand2 routine, CPML-10, B-3 
atandf routine, B-2 
atand routine, CPML-9, B-2 
atanf routine, B-2 
atanhf routine, B-3 
atanh routine, CPML-11, B-3 
atan routine, CPML-9, B-2 


Index 


Base 10 logarithm, CPML-46 
Base 2 logarithm, CPML-46 
Bessel functions, CPML-12 
bessel routine, CPML-12 


Cc 


cabsf routine, B-3 

cabs routine, CPML-13, B-3 

cortf routine, B-3 

cort routine, CPML-14, B-3 

ccosf routine, B-3 

ccos routine, CPML-15, B-3 

cdivf routine, B-4 

cdiv routine, CPML-4, CPML-16, B-4 

ceilf routine, B-4 

Ceiling, CPML-17 

ceil routine, CPML-17, B-4 

cexpf routine, B-4 

cexp routine, CPML-18, B-4 

dogf routine, B-4 

dog routine, CPML-19, B-4 

cmulf routine, B-4 

cmul routine, CPML-20, B-4 

Common logarithm, CPML-46 

Complex absolute value, CPML-13 

Complex data types, 1-2 

Complex division, CPML-16 

Complex exponential, CPML-18 

Complex floating-point data types, 1-3 

Complex functions, 1-3 

Complex multiplication, CPML-20 

Complex natural logarithm, CPML-19 

Complex numbers, 1-3 
absolute value, CPML-13 
cosine of angle, CPML-15 
division, CPML-16 
exponential of, CPML-18 
exponentiation of, CPML-25 
multiplication, CPML-20 
natural logarithm of, CPML-19 
sine of, CPML-26 
square root of, CPML-27 


Index—1 


Complex power, CPML-25 
Complex square root, CPML-27 


Converting floating-point data types, CPML-28 


copysignf routine, B-4 

copysign routine, CPML-21, B-4 

cosdf routine, B-5 

cosd routine, CPML-22, B-5 

cosf routine, B-5 

coshf routine, B-5 

cosh routine, CPML-23, B-5 

Cosine and sine of angle, CPML-57 

Cosine of angle, CPML-22 
hyperbolic, CPML-23, CPML-59 
of a complex number, CPML-15 

cos routine, CPML-22, B-5 

Cotangent of angle, CPML-24 

cotdf routine, B-5 

cotd routine, CPML-24, B-5 

cotf routine, B-5 

cot routine, CPML-24, B-5 

cpowf routine, B-6 

cpow routine, CPML-25, B-6 

csinf routine, B-6 

csin routine, CPML-26, B-6 

csqrtf routine, B-6 

csqrt routine, CPML-27, B-6 

Cube root, CPML-14 

cvt_ftof routine, CPML-28, B-6 


D 


Data types 
conversion of floating-point, CPML-28 
INPUT_ARG TYPE, 2-2 
RETURN_TYPE, 2-2 
types of, 1-2 

Decimal boundary values, A-1 

Denormalized numbers, 1-5 

Domain errors, 1-4 

dremf routine, B-6 

drem routine, CPML-32, B-6 

D_FLOAT data type, conversion to and from, 

CPML-28 


E 


Entry-point names, B-1 
Entry points, 2-2 

erfcf routine, B-7 

erfc routine, CPML-33, B-7 
erff routine, B-6 

erf routine, CPML-33, B-6 
Error functions, CPML-33 
Error handling, 2-3 


Errors 
domain, 1-4 
range, 1-4 


Index—2 


Euclidean distance, CPML-41 
Exceptional arguments, 1-3, 1-4, 2-3 
Exception conditions, 1-4, 2-3 
Exception handler, 1-1 

expf routine, B-7 

expmlf routine, B-7 

expm1 routine, CPML-34, B-7 
Exponent adjustment, CPML-55 
Exponential, CPML-34 

Exponential of a complex number, CPML-18 
exp routine, CPML-34, B-7 


F 


fabsf routine, B-7 

fabs routine, CPML-35, B-7 
finitef routine, B-7 

finite routine, CPML-36, B-7 
Finite value, checks for, CPML-36 
Floating-point complex data types 


for Compag Tru64 UNIX Alpha systems, 1-3 


for OpenVMS Alpha systems, 1-3 
IEEE, 1-3 
types of, 1-3 
VAX, 1-3 
Floating-point data types 
complex, 1-2, 1-3 
conversion to other types, CPML-28 


for Compag Tru64 UNIX Alpha systems, 1-2 


for OpenVMS Alpha systems, 1-2 
IEEE, 1-2 
types of, 1-2 
VAX, 1-2 
Floating-point number conversion 
Big Endian_|EEE_S Floating, CPML-28 
Big Endian_|EEE_T Floating, CPML-28 
Big Endian_|EEE_X_Floating, CPML-28 
CRAY_Floating, CPML-28 
D_Floating, CPML-28 
F_Floating, CPML-28 
G Floating, CPML-28 
H_ Floating, CPML-28 
IBM_Long Floating, CPML-28 
IBM_Short_Floating, CPML-28 
IEEE_S Floating, CPML-28 
IEEE_T Floating, CPML-28 
IEEE_X Floating, CPML-28 
floorf routine, B-8 
floor routine, CPML-37, B-8 
fmodf routine, B-8 
fmod routine, CPML-38, B-8 
fp_classf routine, B-8 
fp_class routine, CPML-39, B-8 
frexpf routine, B-8 
frexp routine, CPML-40, B-8 
F_COMPLEX data type, 1-2 


F_FLOAT datatype, 1-2 
conversion to and from, CPML-28 J 
F_FLOAT_COMPLEX data type, 1-3 


jOf routine, B-9 


jO routine, CPML-12, B-9 
G jlf routine, B-9 
gamma routine, CPML-45 jl routine, CPML-12, B-9 
Generic interface names, 2-2 jnf routine, B-9 
G_FLOAT datatype, 1-2 jn routine, CPML-12, B-9 
conversion to and from, CPML-28 
G_FLOAT_COMPLEX data type, 1-3 L 
H Language-specific routine behavior, 1-1 
Idexpf routine, B-9 
Hexadecimal boundary values, A-1 Idexp routine, CPML-44, B-9 
Hyperbolic arc cosine of angle, CPML-6 Igammaf routine, B-10 
Hyperbolic arc sine of angle, CPML-8 Igamma routine, CPML-45, B-10 
Hyperbolic arc tangent of angle, CPML-11 Igamma_max_float boundary value, A-4 
Hyperbolic cosine of angle, CPML-23 In routine, CPML-46, B-10 
Hyperbolic sine, CPML-58 logl0f routine, B-10 
Hyperbolic sine and cosine, CPML-59 log10 routine, CPML-46, B-10 
Hyperbolic tangent, CPML-62 loglpf routine, B-10 
hypotf routine, B-8 loglp routine, CPML-46, B-10 
hypot routine, CPML-41, B-8 log2f routine, B-10 
H_FLOAT data type, conversion to and from, log2 routine, CPML-46, B-10 
CPML-28 Logarithm 
base 10, CPML-46 
| base 2, CPML-46 
common, CPML-46 
IEEE Big Endian, CPML-28 complex number of a, CPML-19 
IEEE considerations, 1-5 gamma function, CPML-45 
IEEE double-precision complex data type, 1-3 ilogb, compute unbiased exponent, CPML-42 
IEEE doubleprecision datatype, 1-2 Igamma function, CPML-45 
|EEE extended-precision complex data type, 1-3 logb, convert to double-precision floating-point, 
IEEE extended-precision data type, 1-2 CPML-47 
IEEE floating-point data types, conversion to other natural, CPML-46 
data types, CPML-28 logbf routine, B-10 
IEEE floating-point values, identifying the class of, logb routine, CPML-47, B-10 
CPML-39 logf routine, B-10 
IEEE Little Endian, CPML-28 log routine, CPML-46, B-10 
IEEE single-precision complex data type, 1-3 
IEEE single-precision data type, 1-2 M 
ilogbf routine, B-8 
ilogb routine, CPML-42, B-8 Machine numbers, CPML-49 
Input arguments, 2-2 math$acosd f, B-1 
Integer data type, 1-2 math$acosd_g, B-1 
Integers, rounding, CPML-50, CPML-54 math$acosd_s, B-1 
Integral power of 2, converting to, CPML-40 math$acosd_t, B-1 
Invalid argument exception condition, 1-4, 1-5 math$acosd_x, B-1 
Invalid arguments, 1-4 math$acosh_f, B-1 
isnanf routine, B-9 math$acosh_g, B-1 
isnan routine, CPML-43, B-9 math$acosh_s, B-1 
math$acosh_t, B-1 
math$acosh_x, B-1 


math$acos f, B 
math$acos_g, B-1 
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math$acos s, 
math$acos t, 
math$acos_x, 
math$asind_f, 
math$asind_g, B 
math$asind_s, B 
math$asind_t, B 
math$asind_x, B- 
math$asinh_f, B 
math$asinh_g, B 
math$asinh_s, B 
math$asinh_t, B- 
math$asinh_x, B 
math$asin_f, B-2 
math$asin_g, 
math$asin_s, 
math$asin_t, 
math$asin_x, 
math$atan2 f, B-2 
math$atan2_g, B 
math$atan2_s, B- 
math$atan2_t, B 
math$atan2_x, B- 
math$atand2_f, 
math$atand2_g, 
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B-3 
math$atand2_s, B-3 
B-3 


math$atand2 t, 


math$atand2_x, B-3 


math$atand f, B-2 
math$atand_g, B 
math$atand_s, B- 
math$atand _t, B 
math$atand_x, B- 
math$atanh_f, B- 
math$atanh_g, B 
mathgatanh_s, B- 
math$atanh_t, B 
math$atanh_x, B- 
math$atan_f, B- 
math$atan_g, B-2 
math$atan_s, B 
math$gatan_t, B- 
math$atan_x, B-2 


math$cbrt_f, B-3 
math$cbrt_g, B-3 
math$cbrt_s, B-3 
math$cbrt_t, B-3 
math$cbrt_x, B-3 
math$ccos f, B-3 
math$ccos_g, B-3 
math$ccos_s, B-3 
math$ccos t, B-3 
math$ccos_x, B-3 
math$cdiv_f, B-4 
math$cdiv_g, B-4 
math$cdiv_s, B-4 
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math$cdiv_t, B-4 
math$cdiv_x, B-4 
mathé$ceil_f, 
mathé$ceil_g, 
math$ceil_s 
mathé$ceil_t, 
math$ceil_ x, 
math$cexp f, 
math$cexp_g, 
math$cexp_s 
math$cexp t, 
math$cexp_x, 
math$clog_f, 
math$clog_g, 
math$clog_s, 
math$clog_t, 
math$clog_x, 
math$cmul_f, 
math$cmul_g, 
math$cmul_s, 
math$cmul_t, 
math$cmul_x, 
math$copysign _ 
math$copysign_ 
math$copysign__ 
math$copysign _ 
math$copysign _ 
math$cosd_f, 
math$cosd_g, 
math$cosd_s, 
math$cosd_t, 
math$cosd_x, 
math$cosh_f, 
math$cosh_g, 
math$cosh_s, 
math$cosh_t, 
math$cosh_x, 
math$cos_f, B-5 
math$cos_g, B- 
math$cos_s, B- 
math$cos t, B- 
math$cos_x, B- 
math$cotd_f, B-5 
math$cotd_g, B-5 
math$cotd_s, B-5 
math$cotd_t, B-5 
math$cotd_x, B-5 
math$cot_f, B-5 
math$cot_g, 
math$cot_s, B-5 
math$cot_t, B-5 
math$cot_x, B-5 
math$cpow_f, B-6 
math$cpow_g, B- 
math$cpow_s, 
math$cpow_t, 
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math$cpow_x, B-6 
math$csin_f, B 
math$csin_g, B 
math$csin_s, B- 
math$csin_t, B 
math$csin_x, B 


math$csqrt_ g, B 
math$csqrt_s, B 
math$csqrt_t, B 
math$csqrt_x, B 
math$erfc f, B-7 
math$erfc_g, 
math$erfc_s, B-7 
math$erfc t, B-7 
math$erfc_x, B-7 
math$erf_f, 
math$erf_g, 6 
math$erf_s, 
math$erf_t, 
math$erf_x, B-6 
math$expm1 f, B-7 
math$expm1 4g, 
math$expm1_s, B 
math$expm1 t, B 
math$expm1_x, B- 
math$exp f, B-7 
math$exp g, B- 
math$exp s, B- 
math$exp t, B 
math$exp_x, B- 
math$fabs f, B-7 
math$fabs_g, B-7 
math$fabs_s, B-7 
math$fabs t, B-7 
math$fabs_x, B-7 
math$finite f, B-7 
math$finite_g, B 
math$finite_s, B-7 
math$finite t, B 
math$finite_x, B 
math$floor_f, B-8 
math$floor_g, B-8 
math$floor_s, B-8 
math$floor_t, B-8 
math$floor_x, B-8 
math$fp_dass f, B 
math$fp_class_g, 
math$fp_class _s, 
math$fp_class t, 
math$fp_class_x, 
math$frexp f, B- 
math$frexp_g, B 
math$frexp_s, B- 
math$frexp t, B 
math$frexp x, B- 


math$hypot_f, 
math$hypot_g, 
math$hypot_s 
math$hypot_t, 
math$hypot_x, B-3, B- 
math$ilogb f, B-8 
math$ilogb_g, B- 
math$ilogb_s, B- 
math$ilogb t, B- 
math$ilogb_x, B- 
math$jO_f, B-9 
math$jO_g, B 
math$j0O_s, B- 
math$jO_t, B 
math$j0_x, 
math$j1_f, 
math$j1_g, B 
math$jl_s, B- 
math$j1_t, B 
math$jn_f, B- 
math$jn_g, B-9 
math$jn_s, B-9 
math$jn_t, B-9 
math$jn_x, B-9 
math$ldexp f, B 
math$ldexp_g, B-9 
math$ldexp_s, B-9 
math$ldexp_t, B-9 
math$ldexp_x, B-9 
math$lgamma_f, B-10 
math$lgamma_g, B-10 
math$lgamma_s, B-10 
math$lgamma_t, B-10 
math$lgamma_x, B-10 
math$In_f, B-10 
math$In_g, B-10 
math$In_s, B-10 
math$ln_t, B-10 
math$In_x, B-10 
math$logl0_f, B-10 
math$logl10_g, B-10 
math$logl0_s, B-10 
math$log10_t, B-10 
math$logl10_x, B-10 
math$loglp f, B-10 
math$loglp_g, B-10 
math$loglp_s, B-10 
math$loglp t, B-10 
math$loglp_x, B-10 
math$log2_f, B-10 
math$log2_g, B-10 
math$log2_s, B-10 
math$log2_t, B-10 
math$log2_x, B-10 
math$logb_f, B-10 
math$logb_g, B-10 
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math$logb_s, B-10 
math$logb_t, B-10 
math$logb_x, B-10 
math$modf_f, B-11 
math$modf_g, B-11 
math$modf_s, B-11 
math$modf_t, B-11 
math$modf_x, B-11 
math$mod_f, B-8 
math$mod_g, 
math$modc_s, 
math$mod_t, 
math$moc_x, 
math$nextafter_f, B-11 
math$nextafter_g, B-11 
math$nextafter_s, B-11 
math$nextafter_t, B-11 
math$nextafter_x, B-11 
math$nint_f, B-11 
math$nint_g, B-11 
math$nint_s, B-11 
math$nint_t, B-11 
math$nint_x, B-11 
math$pow_ff, B-11 
math$pow_gg, B-11 
math$pow_ss, B-11 
math$pow tt, B-11 
math$pow_xx, B-11 
math$random_|_f, B-11 


B 
B 
B 
B 


math$rem_f, B-6, B-12 
math$rem_g, B-6, B-12 
math$rem_s, B-6, B-12 
math$rem_t, B-6, B-12 
math$rem_x, B-6, B-12 
math$rint_f, B-12 


math$rint_g, B-12 
math$rint_s, B-12 
math$rint_t, B-12 
math$rint_x, B-12 
math$scalb f, B-12 
math$scalb_g, B-12 
math$scalb_s, B-12 
math$scalb t, B-12 
math$scalb x, B-12 
math$sincosd_f, B-12 
math$sincosd_g, B-12 
math$sincosd_s, B-12 
math$sincosd_t, B-12 
math$sincosd_x, B-12 
math$sincos f, B-12 
math$sincos_g, B-12 
math$sincos_s, B-12 
math$sincos t, B-12 
math$sincos_x, B-12 
math$sind_f, B-13 


Index—6 


math$sind_g, B-13 
math$sind_s, B-13 
math$sind_t, B-13 
math$sind_x, B-13 
math$sinhcosh_f, B-13 
math$sinhcosh_g, B-13 
math$sinhcosh_s, B-13 
math$sinhcosh_t, B-13 
math$sinhcosh_x, B-13 
math$sinh_f, B-13 
math$sinh_g, B-13 
math$sinh_s, B-13 
math$sinh_t, B-13 
math$sinh_x, B-13 
math$sin_f, B-12 
math$sin_g, B-12 
math$sin_s, B-12 
math$sin_t, B-12 
math$sin_x, B-12 
math$snan_f, B-9 
math$snan_g, B-9 
math$snan_s, B-9 
math$snan_t, B-9 
math$snan_x, B-9 
math$sqrt_f, B-13 
math$sqrt_g, B-13 
math$sqrt_s, B-13 
math$sqrt_t, B-13 
math$sqrt_x, B-13 
math$tand_f, B-14 
math$tand_g, B-14 
math$tand_s, B-14 
math$tand_t, B-14 
math$tand_x, B-14 
math$tanh_f, B-14 
math$tanh_g, B-14 
math$tanh_s, B-14 
math$tanh_t, B-14 
math$tanh_x, B-14 
math$tan_f, B-13 
math$tan_g, B-13 
math$tan_s, B-13 
math$tan_t, B-13 
math$tan_x, B-13 
math$trunc_f, B-14 
math$trunc_g, B-14 
math$trunc_s, B-14 
math$trunc_t, B-14 
math$trunc_x, B-14 
math$unordered_f, B-14 
math$unordered_g, B-14 
math$unordered_s, B-14 
math$unordered t, B-14 
math$unordered_x, B-14 
math$y0 f, B-14 
math$y0_g, B-14 


math$y0_s, B-14 


math$y0_t, B-14 S 
math$y0_x, B-14 
math$y1 f, B-14 scalbf routine, B-12 
math$y1_g, B-14 scalb routine, CPML-55, B-12 
math$y1_s, B-14 sincosdf routine, B-12 
math$y1 t, B-14 sincosd routine, CPML-57, B-12 
math$y1_x, B-14 sincosf routine, B-12 
math$yn_f, B-15 sincos routine, CPML-57, B-12 
math$yn_g, B-15 sindf routine, B-13 
math$yn_s, B-15 sind routine, CPML-56, B-13 
math$yn_t, B-15 Sine, hyperbolic, CPML-58 
math$yn_x, B-15 Sine and cosine of angle, CPML-57 
max_float boundary value, A-1 Sine of angle, CPML-56 
min_float boundary value, A-1 hyperbolic, CPML-59 
modff routine, B-11 Sine of angle of a complex number, CPML-26 
modf routine, CPML-48, B-11 sinf routine, B-12 
Modulo remainder, CPML-38 sinhcoshf routine, B-13 
sinhcosh routine, CPML-59, B-13 
N sinhf routine, B-13 
sinh routine, CPML-58, B-13 
NaN value, checking for, CPML-43 sin routine, CPML-56, B-12 
Natural logarithm, CPML-46 sqrtf routine, B-13 
complex number of a, CPML-19 sqrt routine, CPML-60, B-13 
nextafterf routine, B-11 Square root, CPML-60 
nextafter routine, CPML-49, B-11 of complex numbers, CPML-27 
nintf routine, B-11 Symbolic constants, 1-4 
nint routine, CPML-50, B-11 S FLOAT datatype, 1-2 
Normalized fractions, converting to, CPML-40 S_FLOAT data type, conversion to and from, 
CPML-28 
O S FLOAT data type, IEEE Big Endian, CPML-28 
S_ FLOAT data type, IEEE Little Endian, 
Overflow exception condition, 1-4, 1-5 CPML-28 
S FLOAT_COMPLEX data type, 1-3 
P 
powf routine, B-11 a 
pow routine, CPML-51, B-11 tandf routine, B-14 
tand routine, CPML-61, B-14 
R tanf routine, B-13 
Tangent of angle, CPML-61 
random routine, CPML-52, B-11 hyperbolic, CPML-62 
Range errors, 1-4 tanhf routine, B-14 
Remainder tanh routine, CPML-62, B-14 
drem function, CPML-32 tan routine, CPML-61, B-13 
modulo, CPML-38 Truncation, CPML-63 
remainder function, CPML-53 truncf routine, B-14 
remainderf routine, B-12 trunc routine, CPML-63, B-14 
remainder routine, CPML-53, B-12 T_FLOAT datatype, 1-2 
return type, 2-2 T_FLOAT data type, conversion to and from, 
Right triangle, hypotenuse of a, CPML-41 CPML-28 
rintf routine, B-12 T_FLOAT data type, IEEE Big Endian, CPML-28 
rint routine, CPML-54, B-12 T_FLOAT data type, IEEE Little Endian, 
Rounding to the nearest integer, CPML-50 CPML-28 
Routine interface, 2-2 T_FLOAT_COMPLEX data type, 1-3 


examples, 2-3 
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Underflow exception condition, 1-4, 1-5 
unorderedf routine, B-14 
unordered routine, CPML-64, B-14 


V 


VAX double-precision complex data type, 1-3 
VAX double-precision data type, 1-2 
VAX single-precision complex data type, 1-3 
VAX single-precision data type, 1-2 
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conformance to, 1-5 
X_FLOAT datatype, 1-2 
X_FLOAT data type, conversion to and from, 
CPML-28 
X_FLOAT data type, IEEE Big Endian, CPML-28 
X_FLOAT data type, IEEE Little Endian, 
CPML-28 
X_FLOAT_COMPLEX data type, 1-3 


Y 


yOf routine, B-14 
yO routine, CPML-12, B-14 
ylf routine, B-14 
yl routine, CPML-12, B-14 
ynf routine, B-15 
yn routine, CPML-12, B-15 


